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ABSTRACT 


The  determination  of  aircraft  flotation  performance  and  operations 
capability  on  semi-  and  unprepared  soil  runways  is  a  critical  factor  in 
developing  forward  area  airfields  in  limited  warfare  situations.  This 
investigation,  which  is  a  part  of  a  continuing  program,  was  directed  at 
defining  landing  gear- soil  interaction  and  developing  the  flotation  criteria 
to  permit  comparative  evaluation  of  the  relative  merits  of  various  landing 
gear  configurations. 

A  basic  aircraft  tire-soil  interaction  equation  relating  the  drag  ratio 
(R/P)  to  sinkage  ratio  (Z/D)  was  developed  for  the  Region  II  velocity 
range  (5  knots  to  40  knots).  The  influence  of  high  velocity  and  multiple 
wheel  configurations  on  flotation  performance  was  determined  on  a 
preliminary  basis.  Empirical  sinkage  prediction  equations  were  developed 
for  predicting  the  sinkage  of  aircraft  type  tires  in  cohesive  and  cohesionless 
soils  with  an  estimated  accuracy  of  ±40%  within  the  90%  confidence  limits. 
The  results  of  the  Single  Wheel  Verification  Tests  are  reported  and  used  to 
verify  the  developed  flotation  analysis  equations. 

An  analytical  approach  to  sinkage  prediction  using  finite  element 
techniques  was  developed  to  give  a  more  rational  approach  to  sinkage 
analysis.  The  soil  was  assumed  to  be  an  elastic-perfectly  plastic  medium. 
The  results  of  this  analytical  approach  as  given  by  the  computer  program 
and  the  results  of  a  test  case  evaluation  are  described  in  detail. 

A  preliminary  Single  Wheel  Relative  Merit  Index  (RMI)  was  developed 
for  permitting  a  comparative  evaluation  of  the  flotation  characteristics 
of  aircraft  tires  on  soil.  The  RMI  was  ubed  to  rate  the  flotation  capacity  of 
aircraft  tires  currently  used  on  cargo,  bomber,  and  fighter  aircraft. 


This  document  is  subject  to  special  export  controls  and  each  transmittal 
to  foreign  governments  or  foreign  nationals  may  be  made  only  with  prior 
approval  of  the  Air  Force  Flight  Dynamics  Laboratory  (FDFM),  Wright- 
Patterson  Air  Force  Base,  Ohio  45433. 
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SECTION  I 
INTRODUCTION 


The  problems  associated  with  the  operation  of  aircraft  on  forward 
area  soil  runways  are  well  documented.  These  problems  include  drag 
resistance  due  to  soil  sinkage,  dynamic  structural  loads  caused  by  sink- 
age  and  ground  roughness  leading  to  structural  damage  and  excessive 
fatigue  accumulations  per  aircraft  mission,  multiple  pass  effects  related 
to  continued  runway  use,  and  limited  runway  lengths  and  widths  for  landing 
and  takeoff  operations. 

In  order  to  develop  the  criteria  for  the  efficient  and  effective 
utilization  of  soil  runways  while  including  adequate  safety  provisions,  it 
is  necessary  to  understand  the  interaction  effects  between  the  tire  and 
soil  during  rolling  and  braking  modes.  In  the  past,  tire-soil  interaction 
studies  have  been  conducted  mainly  in  the  area  of  off-the-road  mobility 
for  military  vehicles.  More  recently,  significant  advances  have  been 
made  toward  defining  landing  gear -soil  interaction  and  developing  the 
flotation  criteria  to  permit  comparative  evaluation  of  the  relative  merits 
of  various  landing  gear  configurations( Flotation  in  this  context 
is  defined  as  the  tire-soil  interaction  and  the  ability  of  the  soil  to  support 
the  tire  load  system. 

The  basic  objective  of  this  program  is  to  analytically  define  landing 
gear -soil  interaction  and  to  develop  a  system  for  rating  the  relative  flo¬ 
tation  capacity  of  landirv  gear  contact  elements  and  landing  gear  systems 
during  aircraft  operations  on  semi-  and  unprepared  soil  runways.  The 
primary  variables  associated  with  evaluating  the  flotation  capacity  of  a 
landing  gear  system  include  drag,  sinkage,  multiple  wheel  effects,  and 
braking  effects.  Secondary  variables  which  influence  flotation  capacity 
include  turning,  landing  impact,  and  point  of  rotation  effects.  The  Rela¬ 
tive  Merit  Index  (RMI)  is  the  system  for  incorporating  these  primary  and 
secondary  flotation  variables  in  a  suitable  manner  to  permit  the  compara¬ 
tive  evaluation  of  the  flotation  capacity  of  landing  gear  types  and  config¬ 
urations.  Phase  l(l)  of  this  program  included  a  survey  of  the  flotation 
problem,  establishment  of  the  critical  flotation  parameters,  and  an  in¬ 
vestigation  of  available  flotation  data  leading  to  the  development  of  a  flo¬ 
tation  analysis  equation. 

Phase  II,  described  in  this  report,  included  the  development  of 
an  empirical  sinkage  prediction  equation,  a  finite  element  based  analytical 
sinkage  prediction  equation,  conducting  the  single-rolling  wheel  verifi¬ 
cation  tests,  refinement  of  the  drag  ratio-sinkage  ratio  relationship,  and 
the  development  of  the  single  wheel  Relative  Merit  Index  (RMI)  system  for 
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defining  comparative  flotation  capacity.  The  Phase  III  portion  of  the  program 
will  include  the  development  and  verification  of  multiple  wheel,  high  velocity, 
and  braking  flotation  effects  and  the  establishment  of  the  finalized  RMi  in¬ 
dexing  system  leading  to  the  development  of  a  final  flotation  criteria  for  de¬ 
fining  aircraft  performance  on  semi-  and  unprepared  soil  airstrips.  A  flow 
chart  indicating  the  overall  flotation  research  program  is  shown  in  Figure  1. 
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Figure  1  Projected  Landing  Gear-Soils  Interaction  and  Flotation  Research  Program 


SECTION  II 

DRAG-SINKAGE  ANALYSIS 


1.  Drag  Ratio-Sinkage  Ratio,  Single  Wheel 

The  drag  forces  exerted  by  the  deformed  ground  on  the  aircraft 
wheel  are  an  important  consideration  in  relation  to  aircraft  thrust  for 
determining  the  takeoff  capability  and  required  runway  lengths  for  air¬ 
craft  operations  on  soil  surfaces.  The  results  of  the  Phase  I  effort 
indicated  that  flotation  performance  of  aircraft  operating  on  soil  runways 
is  directly  related  to  the  sinkage  and  drag  encountered  by  the  aircraft 
tire. 


In  the  Region  II  velocity  range  (5  knots  to  40  knots),  the  drag  has 
been  shown  to  be  directly  proportional  to  the  sinkage^)  as  expressed  by 
the  functional  equation 


R/P  =  f(Z/D) 


(1 


where 

R  =  rolling  drag  force 
P  =  vertical  tire  load 
Z  =  instantaneous  soil  sinkage 
D  =  outside  tire  diameter 

Region  II  has  previously  been  shown  to  be  a  velocity  range  in  which  the 
drag  ratio  (R/P)  and  sinkage  ratio  (Z/D)  remain  relatively  constant(^). 
Subsequent  to  the  development  of  the  initial  drag  ratio-sinkage  ratio 
relationship,  additional  flotation  data  became  available  from  the  Single 
Wheel  Verification  Test  Program  and  the  Lockheed-Langley  Test  Pro¬ 
gram^).  Figure  2,  presents  the  drag  ratio  (R/P)  versus  sinkage  ratio 
(Z/D)  based  on  all  available  flotation  data(^^»  The  flotation 

data  includes  the  full  range  of  soil  type  (cohesive  to  cohesionless)  and 
tire  diameters  from  18"  to  70",  which  encompasses  the  range  of  tire 
sizes  for  military  aircraft.  The  Lockheed-Langley  tests  measured  rut 
depths  rather  than  sinkage  (Z).  Measured  rut  depths  in  the  buckshot 
clay  were  converted  to  instantaneous  sinkage  for  inclusion  in  Figure  2 
using  elastic  rebound  corrections  determined  from  previous  aircraft 
tire  tests  on  buckshot  clay(®)  and  the  results  of  the  Single  Wheel  Veri¬ 
fication  tests  for  which  both  permanent  rut  depth  and  elastic  deformation 
were  measured. 


Drag  Ratio-Sinkage  Ratio,  Single  Wheels  on  Soil 


A  least  squares  linear  fit  of  the  data  given  in  Figure  2,  using 
Equation  (1)  yields  a  drag  ratio-sinkage  ratio  as  given  by 

All  Soil 

=  0.  018  +  3.  23(Z/D)  (2) 

for  0.015^0.  12 

The  positive  R/P  intercept  at  zero  sinkage  is  expected  since  some  rol’lng 
resistance  (friction)  is  present  at  zero  sinkage.  A  closer  evaluation  of 
Figure  2  shows  a  distinct  trend  of  the  sand  drag  ratio  relationship  in 
comparison  to  the  clay  drag  ratio  relationship.  A  least  squares  linear 
fit  of  the  data  for  each  soil  type  gives 

Cohesionless  Soil 

f  =  0.048  +  2.  77(|)  (3) 

Cohesive  Soil  (4) 

§=0.000  +  3.85(|) 

which  indicates  that  aircraft  tires  encounter  more  drag  at  low  sinkages 
when  operating  on  sand  in  comparison  to  clay  but  encounter  less  drag  at 
high  sinkages  on  sand  in  comparison  to  clay.  Based  on  an  analysis  of 
Figure  2,  however.  Equation  (1)  provides  a  good  first  approximation  for 
evaluating  rolling  drag  for  aircraft  tires  on  any  soil  type  w'ithin  the  re¬ 
stricted  range. 

2.  Drag  Ratio-Sinkage  Ratio,  Multiple  Wheel 

The  influence  of  relative  wheel  positions  in  multiple  wheel  config¬ 
urations  must  be  considered  in  evaluating  drag  and  sinkage  effects  for 
aircraft  tires  on  soil.  Two  types  of  multiple  wheel  configurations  are  of 
interest:  dual  (twin)  and  tandem  as  indicated  in  Figure  3.  An  important 
consideration  for  tandem  tires  is  the  extent  to  which  the  rear  tire  follows 
in  the  track  of  the  lead  tire. 

While  dual  tires  with  close  spacing  are,  in  many  instances,  objec¬ 
tionable  for  off-the-road  powered  wheel  vehicles,  close  spacing  of  air¬ 
craft  tires  offer  some  slight  benefits^'  in  reducing  total  rolling  drag 
below  that  which  would  be  encountered  if  each  wheel  operated  as  an  isolated 
wheel.  Although  very  little  twin  wheel  data  is  available.  Figure  4  provides 
an  indication  of  tire  spacing  effects  on  the  resulting  drag  resistance. 
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b  =  TIRE  WIDTH 

n=  FRACTION  OF  TIRE  WIDTH 


TWIN  WHEEL 


TANDEM  WHEEL 


D=TIRE  DIAMETER 

m= FRACTION  OF  TIRE  DIAMETER 


Figure  3  Multiple  Wheel  Landing  Gear  Assemblies 
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Considerable  reduction  in  drag  can  be  expected  in  tandem  tires 
operating  on  soil  when  the  rear  tire  follows  in  the  forward  tire's  track 
since  the  rear  tire  will  undergo  less  sinkage.  Figure  5  shows  this  tandem 
tire  effect  for  clay  soil  as  a  function  of  tandem  tire  spacing.  The  higher 
drag  encountered  at  123  inch  spacing  in  comparison  to  the  60  inch  spacing 
is  probably  related  to  elastic  rebound  of  the  soil  after  passage  of  the  lead 
tire. 


Figures  4  and  5  show  that  a  "configuration  sensitivity"  exists  for 
determining  drag  on  multiple  aircraft  tires  and  indicates  the  need  for  a 
separate  flotation  criteria  for  twin  and  tandem  tire  configurations.  The 
results  of  previous  multiple  tire  studies  conducted  by  WESW  did  not 
recognize  any  difference  between  twin  and  tandem  tires  when  determining 
a  load  adjustment  for  the  equivalent  single  wheel  load. 

The  results,  although  preliminary,  indicate  the  need  for  additional 
research  in  order  to  develop  the  required  flotation  relationships  for  mul¬ 
tiple  tire  configurations. 

3.  High  Velocity 

Aircraft  ground  speeds  range  from  static  conditions  up  to  landing 
and  takeoff  speeds  of  120  knots.  Only  recently  have  these  high  speed  effects 
been  recognized  as  a  separate  phenomena  from  low  speed  performance  in 
analyzing  flotation  effects^'  ).  Figure  6  shows  the  typical  trends  of 
the  drag  ratio  and  sinkage  ratio  versus  horizontal  speed  for  aircraft  tires 
operating  on  soil,  and  indicates  at  least  three  distinct  zones  of  velocity 
response.  At  low  speeds  (less  than  5  knots),  the  drag  ratio  undergoes 
considerable  change  with  velocity  due  to  the  viscous  effects  in  the  soil. 

From  5  knots  to  approximately  40  knots,  the  drag  ratio  remains  relatively 
constant,  indicating  that  the  rate  of  loading  effects  are  no  longer  significant. 
In  the  third  velocity  region,  the  drag  ratio  increases  rapidly,  peaks  at 
what  has  been  termed  the  "planing  speed",  and  then  decreases.  The  "planing 
speed"  is  apparently  a  function  of  the  soil  type,  soil  strength,  and  tire 
inflation  pressure(^). 

There  are  at  least  two  ways  to  consider  the  drag  ratio  response 
in  the  Region  III  velocity  range.  On  the  assumption  that  a  tire  traversing 
soil  at  high  velocities  responds  the  same  as  tires  operating  on  slush, 
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Figure  6  Typical  Drag  Ratio,  Sinkage  Ratio  vs  Horizontal  Velocity 
in  Soil 


Boeing^  utilized  the  drag  equations  developed  for  slush  conditions  and 
interpreted  the  drag  as 

D'  =  KCD][ZV2  (5) 

where  D'  =  additional  drag  due  to  high  speed  effects 

K  =  constant  related  to  soil  density  (p)  and  tire 
width  (b)  =  1/  2pb 

C  =  soil  impingement  drag  coefficient 

V  =  horizontal  velocity 

The  soil  impingement  coefficient,  which  must  be  determined  experimentally, 
has  not  correlated  well  with  the  variables  of  soil  strength  and  tire  inflation 
pressures  in  a  recent  experimental  effort^.  The  drag  defined  by 
Equation  (5)  must  be  added  to  the  drag  developed  in  Region  II  (see  Equation  2) 
to  determine  the  total  drag  acting  on  the  tire  in  the  Region  III  velocity  range. 

An  alternate  approach  would  be  to  study  the  relationship  between  total 
drag  and  instantaneous  sinkage  in  a  manner  similar  to  that  defined  by 
Equation  (1).  If  a  direct  correlation  exists  between  the  drag  and  sinkage, 
then  the  rapid  increase  in  drag  should  be  attributed  more  to  rate  loading 
effects  in  soil  (i.  e. ,  sinkage  response  related  to  rate  of  loading)  than  to 
changes  in  tire  contact  geometry  effects. 

While  only  a  limited  amount  of  high  speed  data  is  available,  Figure  7 
shows  the  results  of  a  recent  experimental  program(^).  For  comparison 
purposes,  both  the  least  squares  fit  of  the  high  velocity  test  data  and  the 
least  squares  relationship  of  Equation  (2)  for  the  Region  II  velocity  range 
is  shown.  The  Region  III  defining  equation  is  given  by, 

|  =  0.038  +  3.  53{|)  (6) 

for  0.01  <  ^  <  0.06 

These  results  indicate  an  approximate  15%  difference  in  the  drag  ratio 
response.  While  sufficient  experimental  evidence  is  not  available  at  this 
time  to  suggest  that  the  drag  ratio-sinkage  ratio  relationship  of  the  form 
of  Equation  (1)  can  be  used  to  interpret  drag  throughout  the  entire  velocity 
range,  it  does  indicate  that  research  efforts  in  studying  the  effect  of  rate 
of  loading  on  soil  sinkage  may  ultimately  explain  the  Region  III  drag  ratio 
response  and  "planing  speed"  effects  rather  than  tire  geometry  change 
effects. 


11 


SINKAGE 


SECTION  in 


SINGLE  WHEEL  VERIFICATION  TESTS 


1.  Purpose 

The  results  of  the  drag  ratio/sinkage  ratio  relationship  shown  in 
Figure  2  includes  the  test  results  from  the  Single  Wheel  Verification  Test 
Program  and  describes  the  basic  rolling  tire-soil  interaction  phenomena. 

This  test  program  was  conducted  to 

1.  provide  additional  confirmation  of  the  preliminary  drag  ratio/ 
sinkage  ratio  relationship  developed  in  Phase  I; 

2.  provide  experimental  data  for  the  development  of  an  aircraft  tire- 
soil  interaction  empirical  sinkage  equation;  and 

3.  provide  experimental  data  for  the  preliminary  verification  of  the 
finite  element  based-analytical  sinkage  prediction  equation  currently  under 
development. 

All  tests  conducted  were  single  pass  only. 

2.  Test  Program 

To  accomplish  these  objectives  a  test  program  was  formulated  as 
outlined  in  Table  I.  The  tires,  loads,  and  soil  parameters  were  selected 
to  provide  a  range  in  sinkage  of  from  1/2  inch  to  3  inches  in  the  sand  and 
clay.  The  following  parameters  were  measured  and  evaluated  in  the 
program: 

wheel  vertical  load 
wheel  drag  load 
tire  deflection 

tire  contact  geometry  (rigid) 

soil  properties  (type,  strength,  and  modulus) 

instantaneous  soil  sinkage 

permanent  rut  depth 

horizontal  wheel  velocity 

3.  T est  Facility 

All  tire-soil  interaction  tests  were  conducted  at  the  U.  S.  Army 
Engineers  Waterways  Experiment  Station  (WES)  Mobility  and  Environmental 
Divisions  Test  Facility.  The  basic  facility  corsists  of  a  row  of  soil  carts 
joined  end  to  end  in  which  processed  soil  is  placed  at  predetermined  conditions 
prior  to  conducting  a  tire  test.  The  test  carriage  dynamometer  is  suspended 
from  an  overhead  rail  system  which  guides  the-  carriage  over  the  row  of 
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Table  1 


Design  Test  Program 


Tire 

Type 


Diameter 

(in) 


Vertical  Soil 
Width  Defl.  Load  Strength 

(in)  (%)  (1^3)  (CD 


Velocity 

(ft/sec) 


Soil 

Type 


7. 00-6-6PR 

ft 

II 

II 

If 

It 

8. 50-10-8PR 

II 

II 

II 

II 

II 

7.  00-6-6PR 

If 

II 

II 

II 

II 

8.  50-10-8PR 

II 

II 

II 

II 


18.5 

t .  8 
»  < 

28 

II 

II 

1 1 

It 

II 

II 

11 

45 

It 

11 

11 

II 

11 

11 

25.  1 

IT) 

. 

oo 

28 

ti 

ti 

1 1 

it 

ii 

11 

ti 

tt 

45 

ii 

ii 

1  1 

1 1 

11 

ii 

18.5 

6.8 

28 

11 

11 

11 

11 

11 

11 

It 

It 

45 

tl 

ti 

11 

It 

11 

11 

25.  1 

8.5 

28 

II 

11 

It 

11 

11 

II 

11 

11 

tt 

11 

11 

11 

II 

It 

II 

It 

11 

II 

II 

11 

45 

II 

11 

11 

11 

11 

11 

If 

If 

It 

It 

11 

11 

11 

II 

11 

300 

500 

700 

500 

1000 

1300 

500 

1000 

1300 

1000 

1500 

2000 

500 

625 

700 

500 

800 

1000 

1000 

1500 

2000 

500 

800 

1000 

1300 

1000 

1500 

2000 

1000 

1300 

1500 


25 

20 

11 

ii 

11 

ii 

It 

it 

11 

ii 

ti 

ii 

25 

20 

11 

n 

11 

n 

11 

1! 

II 

tl 

11 

11 

25 

20 

11 

it 

ti 

tt 

11 

it 

11 

ii 

11 

n 

50 

20 

11 

ti 

II 

1 1 

25 

ii 

II 

II 


50 

II 

It 

25 

II 


Mortar  Sard 

II 

II 

II 

It 

II 

Mortar  Sand 

ii 

it 

it 

ti 

it 

Buckshot  Clay 
n 

ii 

ii 

n 

it 

Buckshot  Clay 

it 

n 

ii 

ii 

it 

n 

ti 

ii 

u 

•i 

ii 

it 
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soil  carts.  The  test  carriage  is  towed  by  a  steel  cable  attached  to  the 
carriage.  Vertical  and  drag  loads  are  measured  using  strain  gage  type 
load  cells  mounted  within  the  dynamometer  system.  Figure  8  shows  a 
view  of  the  overall  test  system  and  Figure  9  shows  a  7.00-6  Type  III 
aircraft  tire  mounted  in  the  carriage  dynamometer.  A  complete  descrip¬ 
tion  of  the  Test  Facility  is  available  in  a  previous  WES  report^). 

4.  Test  Tires 


Two  test  tires  were  selected  for  the  program  which  provided  an 
18  inch  to  25  inch  range  in  tire  diameters.  Figures  9  and  10  show  the 
7.00-6  and  8.  50-10  Type  III  tires,  respectively,  and  Table  2  gives  the 
pertinent  tire  geometry  data.  Each  tire  was  subjected  to  a  series  of 
cyclic  loadings  before  the  tire-soil  interaction  tests  to  eliminate  any 
prestress  effects. 

5.  Soil  Tests  and  Preparation 


The  two  test  soils  selected  for  the  program  were  buckshot  clay 
and  mortar  sand.  The  main  intent  of  the  soil  tests  was  to  provide  an 
accurate  description  of  the  test  bed  soil  and  Its  uniformity  of  placement. 
Soil  tests  conducted  included  moisture  and  density  determinations,  mo¬ 
bility  cone  penetration  resistance  (Cl),  California  Bearing  Ratio  (CBR), 
and  plate  bearing  resistance.  A  complete  description  of  each  soil  test 
is  given  in  Appendix  I. 

The  results  of  the  moisture-density  determinations  are  sum¬ 
marized  in  Table  3. 

In  order  to  establish  a  relationship  between  the  cone  index,  CBR, 
and  plate  bearing  modulus  values  for  use  in  the  sinkage  comparison 
study,  controlled  cone  penetrometer,  CBR,  and  plate  bearing  tests  were 
conducted  for  each  soil  condition.  These  results  as  given  by  the  average 
cone  index  over  0  to  6  inches,  the  average  surface  CBR  at  0.  1  inch  pene¬ 
tration,  and  Bekker  modulus  values  are  summarized  in  Table  4. 

Some  difficulties  were  encountered  in  determining  the  Bekker 
parameters,  since  the  semi-logarithmic  plot  of  the  penetration  resistance 
versus  sinkage  for  each  plate  size  did  not  yield  parallel  lines  as  required 
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Figure  8  Test  Carriage  on  Soil  Carts 


Figure  9  Type  IHf  7.00-6  Test  Tire 


Figure  10  Type  III,  8.  50-10  Test  Tire 
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Table  3 

Moisture-Density  Data 


Test 

Design 

Average  Conditions 

Section 

Soil  Type 

Soil  Strength 

Cl  (0  to  6") 

avg 

Dry  Density 
(pcf) 

Moisture  Content 
(%) 

1 

Buckshot  Clay 

50 

79.9 

36.9 

2 

Buckshot  Clay 

25 

75.  1 

42.  2 

3 

Mortar  Sand 

25 

99.5 

2.5 

Test 

Section  Soil  Type 

1  Buckshot  Clay 

2  Buckshot  Clay 

3  Mortar  Sand 


Table  4 

Test  Bed  Soil  Cor  litions 

Design 

Soil  Strength 

Cl  (0  to  6") 
avg' 

Cl 

avg 

0  to  6" 

25 

25.6 

50 

48.  2 

25 

27.0 

Bekker 

Parameters 


CBR 

kc 

k 

cp 

n 

00 

• 

o 

4.4 

10.6 

0.05 

1.4 

4.0 

22.  2 

0.09 

0.7 

C 

4.0 

0.4 
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by  the  Bekker  sinkage  equation.  Figure  11  shows  the  results  of  a  com¬ 
parison  between  the  average  cone  index  and  the  CBR  based  cn  the  data 
used  in  developing  Table  4.  An  approximate  relationship  between  CIaVg 
and  CBR  as  shown  by  the  dashed  line  in  Figure  1 1  is  given  by 

Cl  (psi)  =  32  CBR  (7) 

avg  r 

which  is  similar  to  the  results  found  from  previous  studies^). 

Undrained  triaxiai  tests  were  also  run  on  undisturbed  samples  of 
the  sand  and  clay  taken  from  the  soil  carts  to  determine  the  cohesion  and 
friction  components  of  the  soil  strength.  These  test  results,  which  are 
described  in  detail  in  Appendix  I,  are  summarized  below  in  Table  5. 

Table  5 


Triaxiai  Test  Results 


Test 

Section 

Soil  Type 

Design 

Soil  Strength 

Cl  (0"  to  6") 

avg 

T r iaxial  Res ults 

c  Cp 

Moisture 

Content 
w,  % 

1 

Buckshot  Clay 

25 

202 

1° 

41.  5 

2 

Buckshot  Clay 

50 

382 

1° 

38.  0 

3 

Mortar  Sand 

25 

0 

36° 

0.  5 

These  results  are  used  in  defining  the  yield  and  plastic  deformation  con¬ 
ditions  in  the  finite  element  sinkage  prediction  equation. 

6.  Test  Results  -  Buckshot  Clay 

A  total  of  14  tire-soil  interaction  tests  were  conducted  on  the  clay 
soil  with  the  8.  50-10  tire  and  6  tests  with  the  7.00-6  tire.  Before  each 
test  run,  cone  penetrometer  tests  were  run  at  two  meter  intervals  to  verify 
the  uniformity  of  soil  conditions.  The  results  of  these  soil  uniformity  tests 
for  each  test  condition  (by  test  number)  are  given  in  Appendix  II.  The  test 
parameters  previously  detailed  were  measured  for  each  test  run.  Drag  and 
vertical  loads  were  measured  continuously  along  the  test  track  by  the 
dynamometer  load  cells  and  recorded  in  oscillographic  traces  which  are 
permanent  records.  Axle  height  was  also  recorded  continuously  along  the 
test  run  and  was  used  for  calculating  the  instantaneous  tire  sinkage(^). 
Instantaneous  soil  sinkage  gages  previously  used  by  WES(®)  were  inserted 
in  a  number  of  test  runs  to  verify  the  sinkage  determinations.  Rut  depths 
were  determined  by  the  difference  between  the  original  soil  surface  and  the 
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final  soil  surface  after  passage  of  the  tire.  The  difference  between  the 
instantaneous  sinkage  and  permanent  rut  depth  is  the  elastic  rebound.  All 
tests  were  conducted  at  an  approximate  horizontal  velocity  of  20  feet  per 
second.  The  complete  results  for  the  tire-soil  interaction  tests  are  sum¬ 
marized  in  Table  6. 

One  way  in  which  increased  flotation  can  be  obtained  is  to  increase 
the  tire  deflection.  Figure  12  shows  this  influence  in  terms  of  the  reduced 
sinkage  and  drag.  Although  applicable  only  to  these  tires  and  soil  conditions, 
a  change  in  tire  deflection  from  28%  to  45%  caused  an  approximated  one- 
third  reduction  in  drag.  Figure  13  shows  the  variation  of  the  drag  ratio  (R/P) 
with  sinkage  ratio  (Z/D)  for  each  tire  and  percent  deflection.  No  clear 
separation  of  the  data  exists  by  tire  deflection.  Equation  (2)  is  also  shown 
on  Figure  13  and  shows  again  the  tendency  toward  higher  drags  at  higher 
sinkages  for  aircraft  tires  on  clay  soil  than  for  the  average  of  sand  and 
clay  (all  soil  line). 

As  indicated  by  Figure  13,  it  is  necessary  to  predict  sinkage  in 
order  to  estimate  drag.  Previous  efforts  (^)  and  the  analysis  in  Section  IV 
of  this  report  have  shown  an  empirical  relationship  between  the  sinkage 
characteristic  (Z/^)  and  the  load-strength  ratio  (a/CIilVg)  where 

l  =  contact  footprint  length  (rigid  surface) 

a  =  ratio  of  vertical  load  (P)  to  rigid  surface  contact  area  (A) 

Cl  =  average  cone  index  over  0"  to  6"  depth 
avg 

Figure  14  shows  the  results  of  the  tire-soil  interaction  data  plotted 
in  this  manner.  It  is  evident  from  these  results  that  as  the  a/CIavg  ratio 
exceeds  approximately  0.6,  the  sinkage  characteristic  increases  more 
rapidly  than  in  the  0.  2  to  0.  6  range.  Since  drag  is  directly  related  to  sink- 
age,  the  drag  will  Increase  in  a  corresponding  manner. 

7.  Test  Results  -  Mortar  Sand 

A  total  of  6  tire-soil  interaction  tests  were  conducted  on  the  sand 
soil  with  the  8.  50-10  tire  and  6  tests  with  the  7.  00-6  tire.  Prior  to  each 
test  run,  the  uniformity  of  soil  placement  was  investigated  by  running  cone 
penetrometer  tests  at  two  meter  intervals  along  the  test  section.  The  re¬ 
sults  of  these  tests  are  given  in  Appendix  II.  As  in  the  clay  tests,  the  basic 
parameters  measured  were  vertical  load,  drag,  sinkage,  and  rut  depth  in 
addition  to  the  tire  and  soil  parameters.  The  complete  results  for  the  tire- 
soil  interaction  tests  on  mortar  sand  are  given  in  Table  7.  All  tests  were 
run  at  a  horizontal  velocity  of  20  feet  per  second. 

The  drag  ratio-sinkage  ratio  results  for  the  sand  are  shown  in 
Figure  15,  together  with  the  least  squares  fit  line  for  all  available  data  as 
defined  by  Equation  (2).  Reference  to  Figures  13  and  15  shows  that  rolling 
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Figure  12  Tire  Deflection  Effects  on  Drag  and  Sinkage, 
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Drag-Sinkage  Data  -  Sand 
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Sand  flowed  back  into  rut 
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aircraff  tires  perform  differently  on  clay  and  sand  soil  in  terms  of  the  drag 
ratio  response  although  these  differences  may  at  most  be  only  ±25%  through 

Z 

the  range  of  0.  02  £  —  <  0.  12. 

Figure  16  shows  the  results  of  the  sinkage  characteristic  versus 
load-strength  comparison  for  sand  and  it  also  illustrates  the  rapid  increase 
in  sinkage  (and  drag)  for  a/CIavg  ratios  in  excess  of  approximately  0.6. 
Similar  reductions  in  drag  for  increasing  tire  deflections  as  shown  for  clay 
in  Figure  12  also  occur  for  aircraft  tires  on  sand. 

8.  Summary 

Whiie  the  Single  Wheel  Verification  Tests  were  limited  in  the  range 
of  test  parameters,  these  results,  as  well  as  other  available  flotation  data, 
substanti  ♦  lhe  basic  tire-soil  interaction  equation  previously  developed^) 
and  subse  /  refined  in  Equation  (2).  These  results  have  also  shown 

that  different  urag  ratio/sinkage  ratio  relationships  exist  for  aircraft  tires 
on  sand  and  clay  soil,  although  Equation  (2)  provides  a  good  approximation 
for  defining  drag  for  all  soils.  The  R/P  versus  Z/D  relationship  shows  that 
if  sinkage  can  be  predicted,  then  drag  can  be  determined  for  aircraft  tires 
operating  on  soil.  Interpretation  of  these  test  results  leading  to  sinkage 
prediction  equations  are  given  in  Sections  IV  and  V. 


SECTION  IV 


SINKAGE  PREDICTION  EQUATIONS 


1.  General  Considerations 

The  results  of  the  Phase  I  and  Phase  II  research  effort  have  shown 
that  flotation  performance  of  aircraft  operating  on  soil  runways  is  directly 
related  to  sinkage  and  drag  encountered  by  the  aircraft  tire.  In  the  Region  II 
velocity  range  as  defined  by  Equation  (2),  and  perhaps  in  the  Region  III 
velocity  range  as  shown  by  Figure  7,  the  drag  has  been  shown  to  be  directly 
proportional  to  the  sinkage  for  rolling  tires.  Based  on  this  relationship, 
the  critical  parameter  in  determining  drag  on  an  aircraft  wheel  is  sinkage. 

If  a  prediction  can  be  made  of  the  anticipated  sinkage,  then  the  developed 
drag  forces  can  be  determined  within  suitable  limits  for  calculating  flo¬ 
tation  performance  of  aircraft  tires  on  soil  runways.  A  previous  reviewH) 
of  presently  used  empirical  sinkage  prediction  equations,  such  as  Bekker^^) 
and  Boeing' ^  which  were  developed  primarily  for  mobility  analysis,  indicated 
that  these  equations  were  not  sufficiently  accurate  to  provide  the  required 
confidence  in  the  prediction  of  flotation  performance.  In  order  to  develop 
more  accurate  sinkage  prediction  techniques  and  to  develop  the  required 
flotation  performance  criteria,  an  empirical  and  an  analytical  (finite  element) 
sinkage  prediction  analysis  were  developed  for  incorporating  those  variables 
and  conditions  specifically  related  to  aircraft  operations  on  soil  runways. 

2.  Sinkage  Prediction  -  Empirical 

Numerous  experimental  and  analytical  efforts^'  have  been 

undertaken  to  develop  load-sinkage  equations  which  adequately  interpret 
soil  properties  as  related  to  soil  deformation  characteristics  under  loading. 
For  rolling  wheels  on  soil,  the  complex  soil  response  involves  elastic, 
viscous,  and  plastic  deformations.  Previously  developed  empirical  sinkage 
prediction  equations  have  mainly  been  directed  at  analyzing  off-the-road 
mobility  of  transport  trucks  and  vehicles  and  armored  vehicles.  While  the 
aircraft  on  soil  problem  shares  many  common  variables  with  the  off-the- 
road  mobility  problem,  sufficient  simplifications  exist  to  permit  the  develop¬ 
ment  of  an  empirically  based  sinkage  prediction  equation  related  primarily 
to  aircraft  type  tire  interaction  with  unprepared  soil  airstrips. 

The  soil  property  measurement  selected  for  referencing  the  soil  is 
the  average  cone  index  over  the  0"  to  6"  depth.  While  most  of  the  tire-sink- 
age  experimental  data  available  used  the  cone  index  as  the  soil  strength 
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measurement,  the  primary  reason  for  selecting  the  cone  index  as  a  ref¬ 
erence  is  based  on  an  analysis  of  the  variables  associated  with  the  cone 

penetrometer.  These  variables  include  P  ,  A  ,  D  ,  % ,  c,  m,  and  Z  where 

c  c  c  c 

P  =  cone  vertical  load  (lbs) 
c 

A  =  cone  base  area  (in3) 
c 

D  =  cone  diameter  (in) 
c 

y  =  unit  weight  of  soil  (lb/in3) 
c  soil  cohesion  (lb/in3) 

Cp  =  soil  angle  of  internal  friction 
=  vertical  penetration  of  cone  (in) 

For  clay  type  soils  (cp=0),  arranging  these  variables  in  a  dimension¬ 
less  array  leads  to 


f 


P/A 


D  ’  y  D 


=  0 


(8) 


Since  the  cone  index  is  normally  defined  at  ZQ  =6"  or  as  an  average  over  a 
6  inch  depth,  Zc/Dc  will  be  constant  and  Equation  (8)  can  be  arranged  as 


P  /A 


(9) 


For  cohesionless  soils  (c=0),  a  similar  analysis  of  variables  leads  to 


f 


(10) 


Again, 


for  Z  /D 
c  c 


being  constant,  the  above  equation  can  be  written  rs 


P  /A 


= 


(ID 


Since  the  cone  penetrometer  has  a  constant  diameter  and  the  cone  index  is 
defined  as  the  ratio  Pc/Ac,  Equations  (9)  and  (11)  indicate  that  cone  pene^ 
tration  readings  taken  in  cohesive  and  cohesionless  soils  will  provide  a 
comparative  evaluation  of  the  soil  (y,c;  or  y  ,  bp) . 
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Some  of  the  drag-sinkage  data  included  in  the  sinkage  analysis  used 
different  soil  parameters  (e.g.,  CBR)  for  soil  reference.  Conversion  of 
these  soil  parameters  to  equivalent  cone  index  were  made  using  techniques 
previously  developed  ^  for  specific  soils  to  provide  a  common  reference 
basis. 


The  primary  concept  behind  the  development  of  the  sinkage  prediction 
equations  was  simplicity  in  the  application  of  the  equations  developed  while 
still  maintaining  reliability  in  the  sinkage  prediction  commensurate  with  the 
known  response  of  soil.  Previous  studies  ^  have  indicated  that  in 

the  Region  II  velocity  range,  the  effects  of  tire  tread,  tire  ply,  soil  viscosity, 
and  soil  inertia  are  not  significant  variables  in  determining  sinkage.  The 
variables  selected  for  inclusion  in  the  analysis  are 


Z ,  i ,  A,  P,  and  Cl 

avg 

where  l  -  tire  contact  footprint  length  (rigid  surface)  at  applicable 
tire  deflection 

Cl  =  average  cone  index  over  0"  to  6"  depth, 
avg 

he  basic  dimensionless  equation  format  selectea  was 


Z 

l 


P/A 

Cl 

avg 


=  f 


a 


Cl 


avg 


where 


Z 

l 


a 


Cl 


=  sinkage  characteristic  (as  distinguished  from  the 
sinkage  ratio,  Z/D) 

-  = load-strength  ratio. 


avg 


(12) 


The  effect  of  drag  on  sinkage  due  to  interaction  effects  was  not  included  in 
the  basic  sinkage  equation  since  the  range  of  interest  for  the  drag  ratio  is 
normally  less  than  0.30  for  effective  aircraft  operation  on  soil  surfaces. 

Cohesive  Soils 


A  limited  amount  of  sinkage  data  was  available  for  aircraft  opera¬ 
tions  on  clay  soil.  The  results  of  the  analysis  ^ which 
include  the  more  recent  Lockheed-Langley  and  Single  Wheel  Verification 
Tests  are  shown  in  Figure  17.  The  data  shown  in  Figure  17  includes  two 
different  clays  (Harpers  Lake  and  buckshot),  tire  diameters  ranging  from 
18  inches  to  70  inches,  tire  deflections  from  15%  to  60%,  and  equivalent 
average  cone  indices  (CBR  converted  to  Cl)  ranging  from  16  psi  to  600  psi 
with  sinkages  from  0.  25  to  4.  0  inches. 
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LOAD-STRENGTH  RATIO,  o/CIavg 

Sinkage  Characteristics  vs  Load-Strength  Ratio,  Cohesive  Soil 


Based  on  an  examination  of  the  data  shown  in  Figure  17,  the 
sinkage  can  be  defined  in  two  different  ranges  as  given  by 


-  =  -0.  03  +  0.  19~ — 

l  Cl 

avg 

for  0.  20  £  - £  0.  60 

avg 


(13) 


and  |  =  -0.11  +  0.33  (14) 

avg 

for  0.60  5  ^ —  <;  0.90 
avg 

The  negative  constant  coefficient  in  Equations  13  and  14  results  from  the 
limited  load-strength  ratio  range  used  in  the  analysis.  Reference  to 
Figure  17  indicates  that  if  the  available  data  is  representative  of  clay  type 
soil  response  under  aircraft  tires,  sinkage  predictions  within  the  90% 
confidence  limits  would  be  within  a  ±40%  accuracy  based  on  the  use  of 
Equations  13  and  14  within  its  restricted  region  of  application. 

Cohesionless  Soils 

(4  14  17) 

Considerably  more  sinkage  data'  ’  *  '  was  available  for  tire 

operations  on  cohesionless  soils.  Unlike  the  clay  data,  however,  much  of 
the  data  from  the  experimental  studies  in  sand  did  not  use  aircraft  type 
tires.  Sufficient  data  was  available,  however,  to  examine  the  effects  of 
such  additional  variables  as  perce  it  tire  deflection  (d)  and  tire  diameter 
to  width  ratio  (D/b)  to  show  that  extensions  of  the  existing  analysis  could 
be  made  to  include  interpretation  of  sinkage  performance  of  aircraft  type 
tires  operating  on  cohesionless  soils. 

To  examine  the  influence  of  tire  deflection,  linear  least  squares 
fits  were  made  for  tire  deflections  of  15%,  25%,  and  35%  for  the  4.00-18 
and  9.  00-14  tires.  These  results  as  shown  in  Figure  18  indicate  no 
consistent  trend  in  the  sinkage-footprint  length  ratio  for  varying  tire 
deflections.  An  examination  was  also  made  of  the  influence  of  the  tire 
diameter-to-width  ratio  on  the  sinkage  characteristic  (Z/f).  Aircraft 
tires  have  a  tire  diameter-to-width  ratio  ranging  from  approximatelv 
2.8  to  3.6.  Figure  19  was  developed  to  show  the  effect  of  the  D/b  ratio 
on  the  operating  level  of  each  tire  as  determined  by  the  variation  of  the 
sinkage  characteristic  (Z  / 1)  with  the  load- strength  ratio.  The  line 
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Figure  18  Sinkage-T ire  Deflection  Analysis,  Cohes 


Figure  19  Sinkage-Tire  Diameter/ Width  Ratio  Analysis,  Cohesionless  Soils 


I 


representing  the  operating  level 
fit  line  for  all  available  data  for 
no  consistent  relationship  of  the 
diameter-to-width  ratio. 


of  each  tire  is  the  linear  least  squares 
that  tire.  Analysis  of  Figure  19  shows 
sinkage  characteristic  (Z  / i)  with  tire 


Based  on  the  results  of  these  comparisons,  available  sinkage 
data  in  cohesionless  soil,  which  includes  the  Lockheed-Langley 

and  Single  Wheel  Verification  tests  presented  herein  were  examined  based 
on  Equation  12.  The  results  of  this  analysis,  as  shown  in  Figure  20, 
includes  three  different  soils,  tire  diameters  ranging  from  18  inches  to 
42  inches,  tire  deflection  from  15%  to  45%,  and  average  cone  indices 
ranging  from  15  psi  to  325  psi  with  corresponding  sinkages  from  0.25  to 
3.0  inches.  The  sinkage  response  is  defined  by  two  ranges  as 


*7 

-  =  -0.  015  +  0.24  - - 

1  “  avg 

for  0.  15  *  —  *  0.70 

avg 


(15) 


and 


|  =  -0.020  +  0.  50  ~ 


avg 


for  0.  70  £ 


a 


Cl 


£  1.00 


avg 


Analysis  of  Figure  20  indicates  that  sinkages  can  be  predicted  based  on  the 
90%  confidence  limits  within  an  approximate  accuracy  of  ±40%  for  the 
indicated  region  of  application. 

In  order  to  provide  an  indication  of  the  range  of  performance  for 
aircraft  tires  operating  on  cohesive  and  cohesionless  soils,  Figure  21  was 
prepared  illustrating  this  range.  It  is  evident  from  reference  to  Figure  21 
that  for  all  other  conditions  being  constant,  sinkages  in  cohesionless  soils 
are  larger  than  in  cohesive  soils. 


3.  Sinkage  Prediction  -  Analytical  (Finite  Element) 


The  empirical  sinkage  prediction  equations,  as  indicated  in  the  pre¬ 
ceding  section,  can  at  best  predict  sinkage  to  within  an  approximate  accuracy 
of  ±40%  based  on  90%  confidence  limits  for  the  restricted  range  of  soil 
properties  used.  Since  using  an  empirical  approach  alone  to  develop  adequate 
sinkage  prediction  equations  for  all  soil  types  and  loading  conditions  is 
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prohibitively  costly  and  time  consuming,  an  analytical  approach  using  the 
finite  element  technique  was  developed  to  give  a  more  rational  approach 
to  sinkage  prediction  and  to  provide  a  long  term  solution  to  sinkage  pre¬ 
diction  for  any  contact  element  in  soil.  The  results  of  this  analytical 
approach,  when  correlated  with  presently  available  experimental  data, 
should  lead  to  more  dependable  sinkage  prediction  for  a  broader  range  of 
soil  properties  and  loading  conditions. 

The  analysis  developed  considers  the  soil  as  an  elastic-plastic 
medium  described  by  the  best  currently  available  plastic  flow  stress- 
strain  relations  and  yield  criterion,  and  the  action  of  the  tire  is  considered 
as  a  pressure  distribution  applied  on  the  surface  of  the  medium.  The 
sinkage  is  then  equivalent  to  the  surface  displacement  due  to  the  applied 
pressure  distribution.  To  obtain  the  surface  displacement,  then,  requires 
the  development  of  the  complete  stress  and  displacement  solutions  of  the 
elastic-plastic  problem. 

The  Vertical  Pulse  Loading  Problem 

Initial  efforts  were  directed  at  determining  the  sinkage  caused  by  a 
rolling  aircraft  wheel  where  the  loading  condition  was  approximated, 
according  to  a  dynamic  analogy  by  the  action  of  a  stationary  wheel 
subjected  to  a  time  dependent  vertical  pulse  load.  The  interface  tire-soil 
contact  stress  distribution  was  approximated  by  a  uniform  pressure  distri¬ 
bution  applied  on  the  surface  of  the  soil  medium  over  a  circular  area 
equivalent  to  the  tire  footprint  area  (see  Figure  22).  The  magnitude  of  the 
pressure  was  taken  to  be  equal  to  the  average  pressure  experienced  by  a 
point  near  the  surface  of  the  soil  when  a  rolling  tire  moves  across  the 
point.  The  dependence  of  the  pressure  at  a  point  in  time  has  been  experi¬ 
mentally  measured  for  land  vehicles  and  has  the  form  of  a  load  pulse 
as  shown  in  Figure  23.  The  peak  pressure  of  the  pulse  is  related  to  the 
vertical  load  on  the  aircraft  tire,  while  the  time  duration  of  the  pulse  is 
related  to  the  aircraft  ground  velocity. 

This  vertical  pulse  is  axisymmetric  and  is  essentially  a  two 
dimensional  problem.  The  conversion  of  the  actual  problem  to  the  vertical 
case  in  using  the  finite  element  method  leads  to  a  simplification  in  the 
solution  over  that  of  the  more  complicated  three  dimensional  rolling  wheel 
problem.  The  three  dimensional  rolling  wheel  problem  will  be  considered 
in  the  subsequent  phase  of  this  research  effoi  • . 
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The  Elastic-Plastic  Soil  Medium 


Recent  works  (19>  20,  2  1)  incjicate  that  soil  media  can  be  simulated 
as  an  elastic-plastic  material.  In  the  present  problem,  the  soil  medium 
was  assumed  to  be  a  homogeneous,  isotropic,  time-independent,  tempera¬ 
ture  independent,  and  elastic-perfectly -plastic  material.  The  elastic 
deformation  was  assumed  to  obey  Hooke's  law,  and  plastic  flow  was  assumed 
to  be  described  by  an  incremental  stre s s- strain  relation  which  is  based  on 
the  normality  flow  rule. 

The  yielding  of  soil  is  a  complicated  phenomenon  and  at  present 
there  is  no  satisfactory  theory  to  account  completely  for  all  the  variables 
in  the  yielding  behavior  of  soil.  The  best  currently  available  criterion  is 
the  Prager-Drucker  Yield  Criterion  22)  which  describes  the  yielding 
in  terms  of  two  parameters,  the  cohesion  and  the  friction  angle  of  the  soil. 
These  two  parameters,  which  are  universally  recognized  soil  parameters, 
can  be  determined  from  triaxial  shear  tests.  While  there  are  other  indica¬ 
tors  of  soil  strength  and  deformation  such  as  the  California  Bearing  Ratio 
(CBR)  and  the  Cone  Index  (Cl),  these  parameters,  which  are  normally  used 
for  empirical  flotation  and  mobility  studies,  cannot  be  interpreted  in  rela¬ 
tion  to  yield  functions  and  plastic  deformation  in  analytical  studies. 

Analytical  Finite  Element  Approach  -  The  Lumped  Parameter 

Iteration  Method 

A  survey  of  past  and  current  literature  indicates  that  both  the  vertical 
loading  and  rolling  wheel  problem  can  be  solved  using  the  finite  element 
approach.  Finite  element  analysis  has  been  used  quite  effectively  in  the 
solutions  of  problems  in  nuclear  blast  effects  and  dynamic  soil  response 
of  strip  footings  These  problems  are  quite  similar  to  the  present 

vertical  pulse  loading  problem.  The  primary  advantages  of  a  finite  element 
type  of  approach  are: 

1.  It  permits  the  inclusion  of  the  plastic  response  of  the  soil. 

2.  It  provides  complete  flexibility  in  applying  surface  geometry 
and  pressure  boundary  conditions  (this  is  an  important  feature  for  extending 
the  work  to  other  contact  elements). 

There  are  two  finite  element  techniques  that  are  applicable  to  the 
vertical  pulse  loading  problem.  They  are  the  direct-stiffness  matrix 
method  ^7)  arKj  the  lumped  parameter  iteration  method  (^6),  The  iatter 
method  was  chosen  for  the  present  problem  because  the  computer  program 
would  be  easier  to  develop  and  fewer  problems  would  be  encountered  with 
computer  capacity,  although  it  would  require  longer  computing  time. 
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In  the  lumped  parameter  iteration  method,  the  continuous  soil 
medium  is  divided  into  lumped  mass  elements  called  mass  points,  and  spring 
elements  called  stress  points  which  can  transmit  both  normal  and  shear 
stresses  between  the  lumped  masses.  These  mass  and  stress  points  are 
interwoven  alternately  as  shown  in  Figure  24.  In  this  model,  the  variables 
are  indicated  only  at  discrete  points.  Consequently,  displacements,  velo¬ 
cities,  and  accelerations  may  be  prescribed  and  solved  for  only  at  the  mass 
points,  while  stresses  and  strains  may  be  prescribed  and  solved  for  only  at 
the  stress  points. 

The  pressure  boundary  conditions  are  applied  on  the  surface,  z  =  0, 
through  a  row  of  fictitious  stress  points  (see  Figure  24).  Outside  the  limits 
of  the  loaded  area,  the  vertical  and  shear  stresses  of  the  fictitious  stress 
points  are  zero.  Within  the  loaded  area,  the  vertical  stress  of  the  fictitious 
stress  points  are  equal  to  the  prescribed  pressure,  and  for  the  present 
problem  the  shear  stresses  are  such  that  the  radial  displacement  on  the 
surface  is  zero.  A  shear  stress  distribution  can  be  prescribed  in  future 
research  efforts  to  take  into  account  effects  such  as  braking  and  turning. 

Using  the  finite  element  model  and  the  method  of  applying  the  boundary 
conditions  described  above,  the  lumped  parameter  iteration  method  obtains 
the  solution  of  the  problem  by  applying  the  prescribed  pressure  according  to 
the  load  pulse  (see  Figure  23)  through  small  increments  of  time,  and  finding 
the  stresses  and  displacements  at  the  discrete  points  of  the  model  for  each 
increment  of  time  by  a  point-by-point  calculation  and  iteration  procedure. 

Governing  Equations 

The  equations  of  continuum  elasticity  and  plasticity,  which  include 
the  equations  of  motion,  quadrature  equations,  yield  and  plastic  equations, 
strain  displacement  relations,  and  the  constitutive  equations,  are  given  in 
this  section  in  the  form  applicable  to  the  finite  element  model.  These 
equations  form  the  basis  for  the  development  of  the  computer  program. 

(Refer  to  the  List  of  Symbols  for  detailed  nomenclature.  ) 

(a)  Dynamic  Equations  of  Motion  -  It  was  shown  by  Ang  (^6)  that  the 
equations  of  motion  for  the  lumped  parameter  finite  element  model  described 
earlier  are  the  same  as  the  differential  equations  of  motion  for  the  corres¬ 
ponding  continuum  model  expressed  in  terms  of  central  finite  difference. 

The  differential  equations  of  motion  for  the  continuum  model  in  the  polar 
coordinates  are: 
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Figure  24  Finite  Element  Model  -  "Lumped  Parameter 
Iteration  Method" 
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These  equations  written  in  central  finite  difference  form  for  a  typical  mass 
point  (i,  j)  of  the  finite  element  model  are:  (The  superscripts  indicate  the 
time  of  the  variables,  and  the  indices  enclosed  in  parenthesis  indicate  the 
spatial  location  of  the  variables  according  to  the  indices  shown  on  Figure  24.  ) 
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For  the  mass  points  that  lie  on  the  axis  of  symmetry,  these  equations  of 
motion  reduce  to: 


ut(i.j)  =  0  (21) 

f  c/d.j+D  - a‘ (i, j-U  2t*  (i+l.j) 

- 2h  - -  +-^h -  (22) 

(b)  Quadrature  Equations  -  The  velocities  and  displacements  at  a 
typical  mass  point  at  time  t  are  approximately  determined  from  the  velocities 
at  time  t-At,  displacements  at  time  t-At,  and  the  accelerations  at  time 
t-At  and  t  by  the  quadrature  equations  which  are  derived  as  follows. 


47 


By  Taylor's  series  expansion, 
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Equations  (23)  through  (26)  can  be  written  as 
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(c)  Prager-Drucker  Yield  Criterion  -  The  criterion  states  that  if 
the  yield  function,  f,  as  defined  below  is  less  than  zero,  the  stress  point  is 
elastic,  and  if  f  is  equal  to  or  greater  than  zero,  the  stress  point  has 
yielded. 


Yield  Function  =  f  =  al  -  k 
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For  the  stress  points  that  lie  on  the  axis  of  symmetry  the  i 
strain-displacement  relations  become,  7’  h 
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For  the  stress  points  that  lie  on  the  first  row  near  the  surface  (that 

is,  j  =2),  relations  for  AeL(i,  j)  and  Ae*(i,  j)  are  the  same  as  Equations  36 

r  ®  t  t 

and  37  or  Equations  42  and  43,  whichever  apply,  but  Ay  (i,  j)  and  A£  (i,j) 

rz  z 

are  defined  differently  as  follows:  (See  Appendix  IV  for  detailed  derivations.  ) 

Incremental  shear-strain  -  For  stress  points  that  do  not  lie  on  the 
axis  of  symmetry. 
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and  for  the  stress  point  that  lies  on  the  axis  of  symmetry. 
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Incremental  vertical  normal  strain  -  If  the  stress  point  is  elastic  or 
unloading  from  plastic  state, 

Ae*(i,  j)  =  ff*(i.  j)  -  o'  At(i,  j)  -j-}—  [a^i,  j)  +  Ae^lJ)]  (47) 


and  if  the  stress  point  is  plastic  and  loading, 
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For  the  stress  points  that  do  not  lie  on  the  axis  of  symmetry, 

d*(i.  J)  =  4  j- 1)  +  q^(i+l ,  j  +  1)  +  at(i- 1,  j- 1)  +  CTfc(i- 1,  j  +  l)l 

Z  *x  Z  Z  Z  Z 

and  for  the  stress  point  that  does  lie  on  the  axis  of  symmetry. 
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and  the  symbols  X,  Q,  B,  and  J  are  defined  by  Equations  55,  64,  63,  and 
33,  respectively. 

(e)  Constitutive  Equations  -  The  incremental  stress  strain  relations 
applicable  to  the  idealized  elastic -plastic  material  for  a  typical  stress  point 
are  shown  below  with  the  (i,j)  indexing  omitted,  since  all  the  variables  are 
referred  to  *he  same  stress  point. 

Elastic  Relations  -  Hooke's  law  written  in  the  incremental  form: 
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Plastic  relations  -  The  plastic  stress-strain  relations  are  based 
on  the  normality  flow  rule  which  states  that  the  plastic  strain  increment 
tensor  is  given  by 


(57) 


where  f  is  the  yield  function  defined  by  Equation  31  and  A  is  a  scalar 
factor  that  depends  on  the  current  stresses  for  perfectly  plastic  material. 
The  final  incremental  stress-strain  relations  for  Prager-Drucker  yield 
criterion  were  derived  by  Christian  (^9)  an(j  are  presented  below  in 
equivalent  forms  (27), 
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0  if  AW  <0,  unloading 

a,  I,  J  are  defined  in  Equations  (32)  through  (34). 


Then,  the  stresses  at  time  t  are: 
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Stress  Correction  for  Perfectly  Plastic  Yielding  -  The  soil 
medium  has  been  assumed  to  be  a  perfectly  plastic  material,  that  is,  there 
is  no  strain  hardening.  This  type  of  material  requires  that  a^ter  yielding 
occurs,  the  yield  function  f  must  be  near  zero  and  not  signif .cantly  above 
zero.  Since  computational  errors  usually  cause  f  to  vary  significantly 
from  zero,  a  form  of  correction  on  the  stresses  are  usually  applied.  In 
the  present  problem,  the  corrections  were  applied  such  that  the  correction 
stress  vector  is  normal  to  the  yield  function  treating  each  component  of  the 
stress  tensor,  that  is,  ar,  Og,  az,  Trz,  and  Tzr  as  independent  variables. 
The  detailed  formulation  of  the  correction  equations  are  presented  in 
Appendix  V  and  the  resulting  equations  are  shown  below. 
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ana  a  ,  a  ,  a  ,  t  are  the  corrected  stresses, 
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Numerical  Procedure  for  the  Development  of  the  Computer 

Program 

For  convenience  in  the  numerical  calculations,  the  above 
governing  equations  were  first  expressed  in  terms  of  dimensionless  variables. 
The  dimensionless  variables  were  formed  in  the  following  manner:  the 
variables  having  a  dimension  of  stress  were  divided  by  the  yield  stress  in 
shear,  k,  (see  Equation  35);  the  variables  having  a  dimension  of  length  were 
divided  by  the  radius,  a,  of  the  applied  surface  pressure  distribution;  and 
time  was  divided  by  the  time  duration  of  the  pulse,  tjj. 

The  pulse  curve  of  the  applied  pressure  on  the  surface  of  the  soil 
was  approximated  by  linear  segments.  For  each  segment,  the  change  in 
pressure  through  an  increment  of  time  was  determined.  The  size  of  the  time 


53 


increment  was  determined  from  criteria  that  will  be  discussed  in  a  subsequent 
section.  For  each  increment  of  time  the  stresses  at  the  fictitious  stress  points 
in  the  plane  (j  =  1)  were  changed  based  on  the  pressure  increments  of  the 
particular  linear  segment  of  the  approximate  pulse  curve.  Then  the  following 
steps  were  performed  starting  with  the  mass  point  at  (i,  j)  =  (3,2): 

(1)  The  accelerations  u  and  w  at  time  t  were  obtained  by  means 
of  the  dynamic  equations  of  motion,  Equations  19  and  20  or  Equations  21  and 
22,  using  the  most  current  stresses  that  were  known  at  time  t  (if  not  known, 
those  at  time  t-At  were  used). 

(2)  The  accelerations  at  time  t  and  the  accelerations,  velocities, 
and  displacements  at  time  t-At  were  then  substituted  into  the  quadrature 
equations,  Equations  27  through  30,  to  give  u,  w,  u,  and  w  at  time  t.  Then 
the  incremental  displacements,  Au*  and  Aw*,  were  obtained  from  Equations  40 
and  41. 

(3)  The  stresses  at  the  three  stress  points  immediately  surround¬ 
ing  the  mass  point  (i,j),  that  is,  the  stress  points  (i,j+l),  i-l,j),  and  i+l,j), 
were  recalculated  according  to  the  following  steps  for  each  stress  point: 

(i)  The  yield  indicating  table  was  checked  to  determine  if 
the  stress  point  had  yielded.  (Initially,  the  table  would  indicate  all  stress 
points  to  be  elastic. ) 

(ii)  The  incremental  strains  at  time  t  were  calculated  using 
the  incremental  strain-displacement  relations.  Equations  36  through  50, 
whichever  set  applied. 

(iii)  The  stress  increments  were  then  calculated  from  the 
incremental  stress  strain  relations.  If  the  stress  point  were  elastic,  the 
elastic  relations.  Equations  51  through  56,  were  used.  If  the  stress  point 
were  plastic,  the  plastic  relations,  Equations  58  through  64,  were  used.  The 
stresses  at  time  t  were  then  calculated  by  Equations  65  through  68. 

(iv)  The  newly  obtained  stresses  were  then  substituted  into 
the  yield  criterion,  Equation  31,  to  check  if  the  stress  point  had  yielded.  The 
result  was  then  recorded  in  the  yield  indicating  table.  If  the  stress  point 
was  plastic,  and  the  yield  function  f  had  exceeded  zero  above  the  limit  at 
which  stress  correction  is  required,  the  newly  calculated  stresses  were 
corrected  by  Equations  69  through  73. 

(4)  Steps  1  through  3  were  repeated  for  the  rest  of  the  mass 
points,  proceeding  from  the  axis  outward  and  then  row  by  row  downward. 

(5)  Using  the  new  stresses  obtained  fer  all  the  stress  points, 
steps  1  through  4  were  repeated  thus  starting  the  iteration  cycle.  This  was 
done  until  the  desired  accuracy  was  reached.  It  was  determined  that  about 
5  or  6  iterations  produced  the  desired  accuracy. 
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(6)  Steps  1  through  5  were  repeated  for  the  subsequent  time 
increments,  in  which  the  applied  pressure  on  the  boundary  was  incremented 
according  to  the  vertical  pulse  load  curve. 

A  flow  diagram  for  the  computer  program  is  shown  in  Figure  25, 
and  a  listing  of  the  computer  program  is  given  in  Appendix  VI. 

Tolerance  for  Yielding,  Unloading,  and  Stress  Correction 

In  the  Prager-Drucker  yield  criterion,  a  stress  point  is  plastic 
when  the  yield  function,  f,  of  that  point  is  greater  than  zero  and  is  elastic 
if  f  is  less  than  zero.  This  means  that  the  yield  surface  has  zero  thickness. 
In  the  numerical  calculation  the  yield  function  of  a  stress  point  fluctuates 
about  zero  during  the  process  of  loading  due  to  computational  accuracy. 
Therefore,  a  tolerance  thickness  of  the  yield  surface  was  used  to  avoid  the 
fluctuation  between  elastic  and  plastic  states  at  a  stress  point.  A  tolerance 
equal  to  1.5%  of  the  yield  stress  in  shear  was  used  above  and  below  zero. 

A  similar  phenomena  occurs  in  the  unloading  from  plastic  state. 
The  incremental  work  done,  AW,  defined  in  Equation  62  also  fluctuates  about 
zero  causing  intermittent  loading  and  unloading  in  the  plastic  state.  To 
avoid  this,  a  tolerance  of  2%  of  the  maximum  incremental  work  done  up  till 
that  time  was  used  for  below  zero  limit. 

Stress  correction  on  the  newly  calculated  stresses  was  applied 
only  if  the  yield  function  is  above  zero  past  a  certain  limit.  In  the  present 
problem,  2%  of  the  yield  stress  in  shear  was  considered  as  the  limit.  It 
should  be  noted  that  all  the  tolerances  used  were  chosen  somewhat 
arbitrarily  but  were  considered  to  be  reasonable  for  this  type  of  problem. 

Selection  of  Mesh  Size  and  Time  Increment 


The  vertical  pulse  loading  is  a  time-dependent  dynamic  problem 
and  consequently,  the  selection  of  the  mesh  size  and  the  selection  of  time 
increment  size  are  interdependent.  There  are  two  criteria  for  selecting  the 
mesh  size  and  the  time  increment.  The  first  criteria  is  related  to  the  con¬ 
vergence  of  the  iteration,  and  the  second  is  related  to  the  stability  of  the 
numerical  calculation. 

The  convergence  criterion  depends  on  the  constant  factor  that 
multiplies  with  the  u’t  and  w* *  terms  in  Equation  28  and  29  of  the  quadrature 
equations.  For  the  quadrature  equation  used,  the  criterion  is  given  by 

Space  Mesh  Size  _  2_h  >  9j_ 

Time  Increment  At  6 
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Figure  25  Flow  Chart  of  Computer  Program 
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where  =  the  velocity  of  elastic  dilatational  wave  propagation  in  the  soil 
medium. 


The  stability  criterion  requires  that 

Space  Mesh  Size  _  2_h  >;> 

Time  Increment  At 

It  was  found  that,  in  order  to  avoid  instability  (^8,  29)  this  ratio  must  be  at 
least  three  times  larger  than  Cj.  In  the  present  problem,  the  space  mesh 
size  was  first  selected  such  that  there  were  at  least  five  fictitious  stress 
points  within  the  radius  of  the  loaded  area.  The  time  increment  was  then 
determined  according  to  the  stability  criterion  given  above.  Because  of 
this  restriction,  the  time  increment  used  was  unavoidably  small  and  conse¬ 
quently  for  load  pulses  that  have  long  time  durations,  the  computer  time 
required  for  the  numerical  calculations  would  be  quite  long. 

Region  of  Solution 

Although  the  problem  to  be  analyzed  is  a  loaded  semi-infinite  half 
space  the  use  of  numerical  procedures  requires  that  the  extent  of  the  region 
influenced  by  the  load  be  restricted  to  a  finite  region.  If  the  boundary  of 
the  mesh  selected  is  not  sufficiently  deep,  an  elastic  rebound  wave  will 
propagate  upward  causing  inaccuracies  in  the  solution  near  the  region  of 
the  loaded  area.  Although  researchers  are  investigating  ways  of  preventing 
this  type  of  reflected  wave  by  modification  of  the  lower  boundary  conditions, 
no  satisfactory  procedure  has  been  developed  to  date.  This  reflected  wave 
condition  was  avoided  in  the  numerical  solution  developed  herein  by  selecting 
the  lower  boundary  sufficiently  deep  to  prevent  any  significant  influence  of 
this  phenomenon  on  the  solution  in  the  region  of  the  loaded  area. 

Test  Case  Evaluation 


A  test  case  having  the  following  typical  soil  and  loading  parameters 
was  analyzed  using  the  finite  element  computer  program  developed. 


Soil  Parameters: 

Density 

Shear  modulus 

Poisson's  ratio 

Cohesion 

F riction  angle 

Yield  stress  in  shear 


Y  =  110  lbs/cu  ft 

G  =  0.  546  x  10®  psf 

V  =  0.4 

c  =  500  psf 
tp  =  20  degrees 
k  =612.3  psf 
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Loading  Parameters: 

Radius  of  the  loaded  area  a  =5.0  inches  f 

Peak  pressure  C7Q  =  10,  000  pif 

Time  duration  of  pulse  t^  =  0.05  sec 

Computational  Parameters: 

Space  Mesh  Size  h  =0.5  inch 

-e 

Time  increment  At  =  25  x  10  sec 

Boundary  size  15.5"  radius  x  15.  5"  deep 

This  test  case  was  run  until  the  applied  pressure  was  at  the  peak  load  and 
a  large  portion  of  the  soil  medium  under  the  loaded  area  had  become  plastic. 
The  plastic  zone  started  along  a  ring  just  under  the  boundary  between  the 
loaded  and  unloaded  area  and  developed  towards  the  axis  along  a  45-degree 
conical  surface  as  the  applied  surface  pressure  increased.  The  development 
of  the  plastic  zone  with  time  and  loading  is  shown  in  Figure  26.  This  plastic 
zone  development  is  consistent  with  other  known  plastic  soil  response,  for 
example,  the  corresponding  plane  strain  static  problem  obtained  by 
Christian  (*9). 


The  vertical  deflections  at  the  surface  of  the  soil  were  plotted 
against  the  radius  from  the  axis  of  the  loaded  area  and  is  shown  in  Figure  27. 

Atj—  =  0.50,  the  applied  pressure  was  sufficiently  large  to  push  some 
td 

portion  of  the  unloaded  surface  above  its  original  level.  At  that  same  applied 
pressure  of  69.5  psi,  the  sinkage  at  r=0  was  comparable  to  the  estimated 
sinkage  based  on  empirical  equations  for  sinkage  prediction.  The  vertical 
deflections  at  the  middle  of  the  applied  pressure  (i.e.  ,  at  r  =  0)  and  at  various 
depths  from  the  soil  surface  were  plotted  against  the  time  of  loading  as 
shown  in  Figure  28. 


Figures  29  and  30  give  the  constant  vertical  stress  contours  at 

—  *  0.  35  and  =  0.  50,  respectively.  The  comparison  of  the  two  figures 

indicates  the  increase  of  vertical  stress  within  the  soil  medium  as  the 
applied  pressure  increases. 


Currently,  the  soil  and  loading  parameters  of  a  few  selected  tests 
from  the  Single  Wheel  Verification  Test  Program  (see  Section  III)  cire  being 
used  in  the  computer  program.  The  selected  tests  were  chosen  such  that 
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DIMENSIONLESS  DEPTH,  z/Q 


l 

I  DIMENSIONLESS  RADIUS,  r/a 


0.2  0.4  0.6  0.8  1.0  1.2  1.4  1.6  1.8  2.0 


Figure  26  Development  of  the  Plastic  Zone  vs  Time  of  Loading 
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Vertical  Deflection  of  the  Surface  vs  Radius  at  Various  Times  of  Loading 


{!•<*)  °*  *3dnss3dd  aanddv 
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DIMENSIONLESS  TIME  t/t 


DIMENSIONLESS  DEPTH,  z/a 


DIMENSIONLESS  RADIUS,  r/o 


|  0.2  0.4  0.6  0.8  1.0  1.2  1.4  1.6 


Figure  29  Constant  Vertical  Stress  Contour  at  Time 
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they  would  be  representative  of  a  range  of  sinkage  in  both  clay  and  sand.  The 
results  of  these  cases  will  be  available  at  a  later  date  in  the  form  of  an 
Interim  Report. 


4.  Comparative  Sinkage  Prediction  Study 


The  results  of  the  Single  Wheel  Verification  Tests  were  used  as  a 
basis  on  which  to  compare  currently  used  sinkage  prediction  equations 
(Bekker  and  Boeing-WES)  and  the  revised  University  of  Dayton  Research 
Institute  (UDRI)  sinkage  prediction  equation.  The  UDRI  sinkage  prediction 
equation  (as  developed  in  Section  IV-2)  is  based  on  laboratory  tire-soil 
test  data  developed  by  WES  and  UDRI,  field  tire-soil  test  data  developed  by 
WES  and  Lockheed-Langley ,  and  aircraft  operational  tire-soil  data  developed 
by  Lockheed  and  Boeing. 


The  Bekker  sinkage  equation  as  given  in  Equation  74  was  used  together 
with  the  soil  properties  data  from  Table  4 


P  = 


(74) 


where  p  =  tire  contact  stress,  psi 

kc,  k^  =  soil  deformation  moduli  (Table  IV) 
n  =  exponent  of  sinkage  (Table  IV) 
b  =  smallest  dimension  of  contact  area,  in. 

Z  =  vertical  sinkage,  in. 

Previously  developed^)  curve  fits  to  the  Boeing  (WES)  sinkage  relationships 
gave  the  sinkage  equations  as 


-  =  0.003  +(CMN) 

for  3£  CMN  £  10 

(75) 

|=  0.  003  +  (SMN)  1,8 

for  5  £  SMN  ^  40 

(76) 

where 

CMN  =  Clay  Mobility  Number  =  ^~~p - (d) 

(77) 

3/8 

SMN  =  Sand  Mobility  Number  =  - — 

(78) 
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and 


Cl  =  average  cone  index  over  0"  to  6",  psi 


G  =  slope  of  cone  index  versus  depth  averaged  over  a  depth 
equal  to  the  tire  width,  psi  per  in. 

b  =  tire  width 

d  =  tire  deflection,  percent 

The  defining  soil  properties  were  developed  from  the  soil  data  as  given  in 
Table  IV.  The  revised  UDRI  sinkage  prediction  equations  are  given  in 
Equations  13  through  16  and  used  soil  property  data  as  given  in  Table  4. 

Figures  31  through  40  show  the  experimental  trend  of  the  sinkage  ratio 
/  D)  versus  increasing  vertical  load  (P)  from  the  Single  Wheel  Verification 
Tests  and  also  the  corresponding  trends  as  predicted  by  the  Boeing  (WES) 
and  UDRI  equations.  The  use  of  the  Bekker  equation,  even  when  accounting 
for  the  increased  tire  contact  area  with  sinkage,  resulted  in  sinkages  which 
in  most  cases  were  an  order  of  magnitude  higher  tha  '  the  measured  sinkage, 
and  consequently,  these  results  are  not  shown  on  Figures  31  through  40. 

The  dashed  line  used  in  the  figures  indicates  the  predictions  are  outside  the 
restricted  range  of  use  of  the  prediction  equations.  Analysis  of  Figures  31 
through  40  indicates  that  the  UDRI  prediction  equations  which  are  based  on  a 
broader  range  of  tire  variables  and  soil  conditions  would  be  more  reliable 
for  predicting  the  sinkage  of  aircraft  type  tires  in  clay  and  sand  types  of 
soil. 

5 .  Summa  ry 

The  results  of  the  sinkage  study  have  shown  that  the  UDRI  sinkage 
prediction  equations  are  presently  the  most  reliable  method  of  predicting  the 
sinkage  of  aircraft  type  tires  in  sand  and  clay  soil  for  the  Region  II  velocity 
range  (5  knots  to  40  knots).  These  empirically  developed  equations  can  be 
used  to  predict  sinkages  with  an  approximate  accuracy  of  ±40%  with  90% 
confidence  within  the  restricted  range  of  analysis  (see  Equations  13  through 
16).  This  restricted  range  of  application  would,  however,  cover  most 
conditions  encountered  by  aircraft  operating  on  soil.  The  effects  of  high 
velocity  (Region  III)  on  sinkage  are  not  adequately  defined  at  this  time  and 
considerable  emphasis  in  research  efforts  should  be  placed  in  studying  high 
velocity  effects  on  tire-soil  interaction. 

The  finite  element  approach  to  sinkage  analysis  using  elastic-plastic 
soil  response  should  ultimately  provide  a  reliable  method  for  extending 
tire-soil  interaction  studies  to  other  soil  types  and  other  contact  elements 
without  undertaking  extensive  experimental  efforts.  The  results  of  the  com¬ 
parisons  between  analytical  (finite  element)  sinkages  and  those  determined  in 
the  Single  Wheel  Verification  Tests  which  will  be  available  in  a  subsequent 
report  should  provide  considerable  insight  into  tire-soil  interaction  phenomena. 
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VERTICAL  TIRE  LOAD,  P(lbs) 


Figure  32  Sinkage  Ratio  vs.  Vertical  Tire  Load,  8.  50-10  Tire 
on  High  Strength  Clay,  d  -  45% 
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Figure  33  Sinkage  Ratio  vs.  Vertical  Tire  Load,  8.  50-10  Tire 
on  Low  Strength  Clay,  d  =  28% 
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Figure  34  Sinkage  Ratio  vs.  Vertical  Tire  Load,  8.  50-10  Tire 
on  Low  Strength  Clay,  d  =  45% 
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Figure  35  Sinkage  Ratio  vs.  Vertical  Tire  Load,  7.00-6  Tire  on 
Low  Strength  Clay,  d  =  28% 
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Figure  36  Sinkage  Ratio  vs.  Vertical  Tire  Load,  7.00-6  Tire  on 
Low  Strength  Clay,  d  =  45% 
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Figure  39  Sinkage  Ratio  vs.  Vertical  Tire  Load,  7.00-6  Tire 
on  Sand,  d  =  28% 


VERTICAL  TIRE  LOAD,  P(lbs) 

Figure  40  Sinkage  Ratio  vs.  Vertical  Tire  Load,  7.00-6  Tire 
on  Sand,  d  =  45% 
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SECTION  V 


SINGLE  WHEEL  RELATIVE  MERIT  INDEX  (RMI) 


1.  Preliminary  Development 

The  primary  variables  associated  with  evaluating  the  flotation  capacity 
of  a  landing  gear  system  includes  drag,  sinkage,  multiple  wheel  effects,  and 
braking  effects.  Secondary  variables  which  influence  flotation  capacity 
include  turning,  landing  impact,  and  point  of  rotation  effects.  The  Relative 
Merit  Index  (RMI)  is  a  system  for  incorporating  these  primary  and  secondary 
flotation  variables  in  a  suitable  manner  to  permit  the  comparative  evaluation 
of  the  flotation  capacity  of  landing  gear  types  and  configurations. 

The  results  of  the  preliminary  development  of  the  basic  tire-soil 
interaction  equation  (1)  and  empirical  sinkage  prediction  equation  (2)  were 
used  to  develop  the  flotation  criteria  for  rating  the  performance  of  aircraft 
tires  operating  on  soil.  This  preliminary  RMI  has  been  developed  in  a 
separate  report  form  and  is  currently  under  review  by  governmental  agencies. 
Copies  of  this  report  are  available  either  from  AFFDL/FDFM  or  the  University 
of  Dayton.  The  preliminary  drag-sinkage  interaction  effects  were  defined  by 


R/P  =  0.015  +  3.09  (|)  (79) 

for  the  Region  II  velocity  range,  and  the  preliminary  sinkage  prediction 
equations  which  were  developed  for  sand  and  clay  are  given  by 


Cohesive  Soil 


Z 

i 


-0.03  +  0.19^j— 
avg 


for  0.20  ■— - s  1.  05 

avg 

where  l  =  contact  footprint  length  (rigid  surface)  at  the 
applicable  tire  deflection 

CIavg  =  average  cone  index  over  0"  to  6"  depth 

a  =  P/A 


(80) 


and 


A  =  tire  contact  area  (rigid  surface)  at  the  applicable 
tire  deflection 
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Cohesionless  Soil 


-  =  -0.01  +  0.23  (81) 

1  avg 

for  0.  015  ^  --g  -  -  <:  1.20 
avg 

Combining  Equations  79>  80,  and  81  leads  to  the  single  wheel  drag  ratio 
relationship  as  given  by 

Cohesive  Soil 

R/P(single  wheel)  =0.  015  +  (0.587^ - -0.093)-^  (82) 

avg 

for  0.  20  <;  ~ —  <;  1.20 
avg 

Cohesionless  Soil 

R/P(single  wheel)  =0.015  +  (0.711^ — -0.031)^  (83) 

avg 

for  0.  15  *  ~ ^  1-  05 
avg 

These  basic  tire-soil  interaction  equations  and  sinkage  prediction  equations 
have  subsequently  been  revised  as  detailed  in  Sections  II  and  IV  based  on 
more  recently  available  flotation  data. 

Equations  82  and  83  were  combined  with  aircraft  tire  geometry 
variables  to  develop  computer  programs  for  evaluating  the  single  wheel  drag 
ratio  relationship.  These  programs,  which  are  given  in  Appendix  III,  form 
the  basis  of  the  preliminary  rolling  single  wheel  comparative  flotation 
evaluation. 

2.  Comparative  Flotation  Evaluation 

The  RMI  system  at  present  has  been  subdivided  into  three  basic 
categories  to  define  the  flotation  capacity  of  aircraft  tires  when  operating  on 
soil  runways.  These  categories  are: 

a.  Takeoff  operations 

b.  Landing  operations 

c.  Runway  deterioration 
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Aircraft  tires  currently  in  use  on  such  aircraft  as  C-119,  C-123,  C-135, 
C-141,  C-5A  and  others  were  rated  using  the  flotation  capability  computer 
program  and  the  defining  index  for  the  tires  operating  on  a  high  strength 
and  a  low  strength  cohesive  and  cohesionless  soil. 

The  preliminary  index  system  associated  with  each  of  these  categories 
are  defined  as: 

Flotation  Index  (FI) 

The  Flotation  Index  is  designed  to  rate  comparatively  the  capability 
of  landing  systems  in  minimizing  the  drag  load  per  unit  vertical  load  for 
specified  operating  conditions  in  takeoff  operations.  The  defining  FI 
equation  is 


FI  =  N[R/P(single  wheel)]  (84) 

where  N  =  the  number  of  contact  elements  in  landing  gear  system.  The 
results  of  the  Flotation  Index  ratings  are  shown  in  Figures  41  through  46. 

Drag  Index  (DI) 

The  Drag  Index  is  designed  to  rate  comparatively  the  capability  of 
landing  gear  systems  in  maximizing  the  drag  load  for  specified  operating 
conditions  in  landing  operations.  An  additional  variable  present  in  the 
landing  mode  is  the  percent  slip  (S)  of  the  braked  tire.  This  braking  con¬ 
dition  leads  to  a  greater  drag  resistance  and  sinkage  for  the  aircraft  tire. 
The  two  basic  parameters  defining  braked  tire-soil  interaction  are: 

—  (S)  =  ratio  of  the  braking  drag  force  to  the  vertical  load  as 

a  function  of  slip;  and 

Z 

—  (S)  =  ratio  of  instantaneous  tire  sinkage  to  the  vertical  load 

as  a  function  of  slip. 

Combining  these  two  parameters  gives  the  defining  Drag  Index  relationship 
as 


DI 


V 


P  ^  rb( 

- -  =  —  (S);  single  wheel 

^(S)  z 


(85) 
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SINGLE  WHEEL  FLOTATION  INDEX  (FI) 

LOW  STRENGTH  SOIL-BOMBER  TIRES 

R-DRAG  FORCE  (#) 


Figure  42 


Flotation  Index,  Aircraft  Bomber 


Tires  on  Low  Strength  Soil 
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SINGLE  WHEEL  FLOTATION  INDEX  (FI) 
LOW  STRENGTH  SOIL-FIGHTER  TIRES 


Flotation  Index,  Aircraft  Fighter  Tires  on  Low  Strength  Soil 


SINGLE  WHEEL  FLOTATION  INDEX  (FI) 

HIGH  STRENGTH  SOIL- CARGO  TIRES 


44  Flotation  Index,  Aircraft  Cargo  Tires  on  High  Strength  Soil 


SINGLE  WHEEL  FLOTATION  INDEX  (FI) 
HIGH  STRENGTH  SOIL-FIGHTER  TIRES 


otation  Index,  Aircraft  Fighter  Tires  on  High  Strength  Soil 


The  Drag  Index  is  then  designed  to  rate  comparatively  the  capability 
of  landing  gear  systems  in  maximizing  the  braking  drag  force  per  unit  instan¬ 
taneous  vertical  sinkage  for  specified  operating  conditions  in  landing  operations. 
Since  an  adequate  braked  tire-soil  interaction  theory  substantiated  by  experi¬ 
mental  verifications  has  not  been  developed  at  the  present  time,  no  Drag  Index 
ratings  were  made  on  existing  aircraft  tires. 

Operations  Index  (OI) 

The  Operations  Index  is  designed  to  rate  comparatively  the  destructive 
(deterioration)  effects  on  soil  runways  of  various  landing  gear  systems  as 
determined  by  the  sinkage  per  unit  vertical  load  for  specified  operating  con¬ 
ditions  in  both  landing  and  takeoff  operations.  The  defining  Ol  equation  is 

OI  =  Z/P  (first  pass)  (86) 

where  Z/P  =  sinkage  factor 

The  results  of  the  Operations  Index  ratings  are  shown  in  Figures  47  through 
52. 


In  addition  to  the  computer  program  for  determining  flotation  capacity, 
a  nomograph  was  developed  for  evaluating  the  Flotation  Index  (FI)  and 
Operations  Index  (OI).  This  nomograph  is  shown  in  Figure  53,  and  an 
example  illustrating  its  use  is  given  below.  Working  Bize  copies  of  the 
nomograph  are  available  from  either  AFFDL/FDFM  or  the  University  of 
Dayton. 


Example  of  Use  of  Nomograph 

Compare  the  relative  flotation  characteristics,  as  defined  by  the 
Flotation  Index  and  the  Operations  Index,  of  the  56  x  16  (Type  VII)  tire  to 
the  40  x  14  (Type  VII)  tire  for  both  tires  operating  at  a  nominal  tire  deflec¬ 
tion  (d)  of  32%. 


Known  Information 


-  Cohesive  Soil,  Cl  =  300  psi 

avg 

-  the  56  x  16  tire  has  a  known  bearing  area  of  340  in3 
at  d  =  40%  and  P  (nominal)  =  76,  000  lbs. 

-  the  40  x  14  tire  has  a  known  bearing  area  of  97  in2 
at  d  =  20%  and  P  (nominal)  =  30,  500  lbs. 
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SINGLE  WHEEL  OPERATION  INDEX  (01) 

LOW  STRENGTH  SOIL-CARGO  TIRES 
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Figure  47 


Operations  Index, 


Aircraft  Cargo  Tires  on  Low  Strength  Soil 
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SINGLE  WHEEL  OPERATION  INDEX  (01) 

LOW  STRENGTH  S01L-B0M8ER  TIRES 


[ 

COHESIVE 

1 

|  NOMINAL  TIRE  DEEL. 

+  IO%TIREDEFL.  | 

49X17(26)- 

56X16(32)- 

36X11(22)- 

24X5  502)  - 

-15X16(24) 

-44X13(26) 

-32X88(22) 

-26X6.604) 

7 

49X17(26)— 

56X16(32)- 

36X11(22)- 

26X66(14) -h 

-15X16(241 

-44X13(261 

-32X8  8(221 

-24X5.5U2) 

7 

Z*VERTICAL  SINKAGE  (IN) 
P»VERTICALTIRE  LOAO(#) 


;  COHESIONLESS 

|  NOMINAL  TIRE  DEFL. 

+10%  TIRE  DEFL. 

49X17(26)- 

-49X17(26) 

56X16(32)- 

-15X16(24) 

56X16(24)— 1 

-56X16(32) 

-44X13(26) 

-44X13(261 

36X0(22  - 

36X11(22)- 

-32X8.8(22) 

32X88(22)- 

-26X6  604) 

Z 

26X6604)- 

7 

-24X5  502) 

7 

-24X55(12) 

Figure  48  Operations  Index,  Aircraft  Bomber  Tires  on  Low  Strength  Soil 


SINGLE  WHEEL  OPERATION  INDEX  (01) 

LOW  STRCNGTH  SOIL-FIGHTER  TIRES 


01  * 

Z/P 

COHESIVE 

L* '0~\ 

NOMINAL  TIRE  OEFL. 

[  +10%  TIRE  DEFL.  ] 

1.0  - 

- 

-36X11(22) 

II  - 

1.2  - 

- 

36X11(221- 

1.3  - 

'  A 

1.4  - 

24X77(14) — 

-30X88(22) 

1.5  - 

1.6  - 

1.7  - 

1.8  - 

1.9  - 

2  0  - 

'increasing' utilization 

24X77(14)- 

25X67504)- 

-30X8  8(221 

25X6  75(141- 

-18X55(12) 

2  1  - 

- 

24X55821- 

2.2  - 

- 

2.3  - 

- 

-18X5502) 

2.4  - 

- 

24X55(121- 

-18X44(8) 

2.5  - 

- 

z 

72.5- 

7 

/ 

18X44(61- 

7 

x' 
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Z 'VERTICAL  SINKAGE  (IN) 
P'VERTICALTIRE  LOADU) 


|  COHESIONLESS 

jNOMINAL  TIRE  DEFL. 

+10%  TIRE  DEFL 

36X111221- 

36XIK221- 

-30X88(22) 

24X7704)- 

-30X8  8(22) 

24X7704)- 

-25X6.7504) 

Z 

18X55(12}— 

24X55(121- 

7 

-25  X  6.7504) 

-18X44(8) 

/. 

18X55(121- 

7 

-24X5  5(12) 

-18X4  4(81 

Figure  49  Operations  Index,  Aircraft  Fighter  Tires  on  Low  Strength  Soil 
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SINGLE  WHEEL  OPERATION  INDEX  (01) 

HIGH  STRENGTH  SOIL-FIGHTER  TIRES 


52  Operations  Index,  Aircraft  Fighter  Tires  on  High  Strength  Soil 


Solution  Nomograph 


-  56  x  16  tire,  D  =  55.7" 

(1)  dj/dg  =  40/32  =  1.  25 

(2)  Intersect  340  ina(A1)  with  d1/ds  =  1.  25  which  gives 


A2  =  280  in2 


(3)  Determine  ^r 


avg 


76,  000 
300 


=  253 


(4)  Intersect  280  in2(A?)  with  — — 

v-»I 


=  253 


avg 


which  gives 


a 


Cl 


=  0.  90 


avg 


(5)  Move  horizontally  to  intersect  the  "cohesive"  soil  line 

2 

and  down  vertically  to  get  —  =  14.  5  x  10"2 

(6)  Compute  l  for  Type  VII  by 

«  =  0.36DVd(4.76-d)  =  0.36x44. 7 Vo. 32(4^7T-3732)  =  23.  8". 

z 

(7)  Intersect  —  =  14.5x  10-2  with  l  -  23.  8"  to  get  Z  =  3.  5". 

Xj 

(8)  Intersect  Z  =  3.  5"  with  D  =  55.  7"  which  gives 


FT,  ,,  =  R/P  =  0.  20 
56x16 

,  3.  5" 

OI56x  16  ~  Z/P  =  76,  000# 


=  4.  6  x  10 


-s 


-40  x  14  tire,  D  =  39.  3" 

(1)  dj/dg  =  20/32  -  0.63 

(2)  Intersect  97  ins(A1)  with  dx  /ds  -  0.63  which  gives 
As  =  155  in2 

(3)  Determine  — 2 -  IP' 

Cl  “  300 

avg 

(4)  Intersect  155  ina(A,,)  with  — ~ —  =  102 

I 

avg 

which  gives  — — ■ - —  =  0.69 
O  I 

avg 
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(5)  Move  horizontally  to  intersect  the  "cohesive  soil"  line 

2 

and  down  vertically  to  get  —  =  10.  5x  10“2 

(6)  Compute  l  for  Type  VII  by 


t.  =  0.36D)/d(4.76-d)  =  0.36  x  39.T\/o. 32(4.76-0. 32)  =  16.  9' 

Z 

(7)  Intersect  —  =  I0.5xl0-2  with  (,  =  16.9"  to  get  Z  =  1.  7" 

(8)  Intersect  Z  =  1.  7"  with  D  =  39.  3"  which  gives 

FI  =  R/P  =  0.15 

40  x  14 


oi4o*m=z/p 


These  results  indicate  that  the  40  x  14  tire  which  would  encounter  25%  less 
drag  per  unit  vertical  load  than  the  56  x  16  tire  has  a  higher  flotation  capacity. 
The  56  x  16  tire  would,  however,  cause  less  runway  deterioration  per  unit 
vertical  load  than  the  40  x  14  tire  as  shown  by  the  OI  ratings. 

3.  Summary 

The  results  of  aircraft  tire-  soil  iteraction  studies  have  led  to  the 
development  of  the  Single  Wbeel-RMI  which  is  a  system  for  evaluating  the 
relative  flotation  capacity  of  aircraft  tires  when  operating  on  soil  runways. 

The  RMI  has  been  presented  in  both  Nomographic  and  Rating  Chart  forms. 

The  potential  uses  of  such  charts  by  landing  gear  engineers  include: 

a.  Design  of  landing  systems  and  contact  element  systems  for 
reduced  sinkage  and/or  drag. 

b.  Comparative  tire  flotation  studies. 

c.  Comparison  of  the  performance  of  aircraft  tires  on  different 
soil  types. 

d.  Determination  of  increased  flotation  with  increasing  tire 
deflection. 

The  preliminary  single  wheel  criteria  and  index  system  developed  to 
date  does  not  include  the  effects  of  multiple  tires  (twin  and  tandem)  and 
braking.  Additionally,  the  influence  of  ground  roughness  and  multiple  passes 
will  ultimately  be  included  such  that  the  RMI  system  can  be  utilized  not  only 
for  comparative  flotation  evaluation  but  also  as  an  absolute  determination  of 
the  suitability  of  a  site  for  aircraft  operations.  As  current  research  effort 
provides  the  required  information  for  evaluating  each  of  these  effects,  the 
RMI  will  be  modified  to  reflect  these  changes. 
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SECTION  VI 


CONCLUSIONS  AND  RECOMMENDATIONS  FOR  RESEARCH 


1.  Conclusions 

The  results  of  the  landing  gear-soil  interaction  and  flotation  criteria 
research  effort  have  shown  that: 

(1)  There  are  at  least  three  velocity  regions  for  which  the  drag  ratio 
exhibits  distinct  characteristics. 

(2)  In  the  Region  II  velocity  range  (*>  knots  to  40  knots),  the  basic 
tire-soil  interaction  as  defined  by  the  drag  ratio  (R/P)  and  sinkage  ratio 
(Z/D)  is  given  by  Equation  2  for  all  soils  although  there  ?re  some  differences 
in  the  performance  of  aircraft  tires  on  sand  and  clay  as  shown  by  Equations 

3  and  4. 

(3)  In  the  Region  HI  velocity  range  (40  knots  to  approximately  70 

to  90  knots),  the  relationship  between  drag  ratio  (R/P)  and  sinkage  ratio  (Z/D) 
has  not  been  adequately  defined;  however,  either  Equation  5  or  6  can  be  used 
on  a  preliminary  basis. 

(4)  Although  considerable  work  remains  to  be  done  on  defining  multiple 
tire  flotation  performance,  preliminary  results  indicate  (a)  that  little  reduc¬ 
tion  in  the  drag  ratio  (R/P)  per  lire  will  occur  for  twin  wheel  configurations, 
and  (b)  that  considerable  reduction  in  drag  ratio  (R/P)  per  tire  can  be  expected 
for  tandem  wheel  configurations  where  the  rear  tire  follows  in  the  track  of 

the  lead  tire. 

(5)  The  empirical  sinkage  prediction  equation  developed  in  this  report 
can  be  used  for  estimating  aircraft  tire  sinkages  in  sand  and  clay  type  soils 
with  an  approximate  accuracy  of  ±40%  within  the  90%  confidence  limits. 

(6)  The  finite  element  based  sinkage  prediction  analysis  described  in 
Section  IV  shows  promise  of  providing  a  rational  (rather  than  empirical) 
method  of  sinkage  prediction  which  could  ultimately  provide  a  better  insight 
into  the  tire-soil  interaction  phenomena  and  lead  to  more  accurate  sinkage 
analysis. 

The  basic  tire-soil  interaction  equation  (Equation  2),  together  with 
the  sinkage  prediction  equations  (Equations  13  through  16)  were  used  to 
develop  the  Single  Wheel  Relative  Merit  Index  (RMI)  which  is  currently  under 
review  by  cognizant  government  agencies.  This  is  the  first  step  leading  to 
the  RMI  of  landing  gear  systems. 
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2. 


Recommendations  for  Research 


The  main  thrust  of  the  immediate  future  landing  gear-soil  interaction 
and  flotation  criteria  research  effort  should  be  directed  at  the  following 
research  deficient  areas. 

(1)  High  speed  (Region  III)  tire-soil  interaction  including  drag  and 
sinkage. 

(2)  Multiple  tire  "configuration  sensitivity"  related  to  drag  and 
flotation  criteria. 

(3)  Flotation  evaluation  of  the  influence  of  braking,  turning,  and 
impact  on  aircraft  performance  on  soil. 

(4)  Additional  evaluation  of  the  failure  criteria  for  defining  the 
capability  of  aircraft  to  operate  on  soil  runways. 

Long  range  flotation  research  effort  should  include  studies  of: 

(1)  Multiple  pass  (operation)  effects  as  related  to  runway  deterioration 
and  capability  of  aircraft  to  operate. 

(2)  Influence  of  runway  roughness  on  the  flotation  performance  of 
aircraft  on  soil. 


9.1 


REFERENCES 


1.  Kraft,  David  C.  ,  Analytical  Landing  Gear-Soils  Interaction,  Phase  I, 

AFFDL-TR-68-88,  Air  Force  Flight  Dynamics  Laboratory,  Wright- 
Patterson  AFB,  Ohio,  August  1968. 

2.  Kraft,  David  C.  and  Luming,  Henry,  Analytical  Landing  Gear-Soils 

Interaction,  Phase  II,  Interim  Report,  University  of  Dayton  Research 
Institute,  Dayton,  Ohio,  January  1969. 

3.  Richmond,  L.  D.  ,  DeBord,  K.  G.  ,  and  Fuller,  J.  R.,  Aircraft  Dynamic 

Loads  from  Substandard  Landing  Sites,  Phase  I,  Interim  Technical 
Report  D6-16190,  The  Boeing  Company,  Renton,  Washington, 
November  1965. 

4.  Crenshaw,  B.M.  ,  Lockheed  Georgia  Co.,  Marietta,  Georgia,  test  data 

furnished  to  David  C.  Kraft,  University  of  Dayton,  May  1969,  Subject: 
Aircraft  Landing  Gear  Dynamic  Loads  from  Operation  on  Clay  and 
Sandy  Soil. 

5.  Frietag,  Dean  R.  ,  Wheels  on  Soft  Soils,  an  Analysis  of  Existing  Data, 

Technical  Report  No.  3-670,  U.  S.  Army  Engineer  Waterways  Experi¬ 
ment  Station,  CE,  Vicksburg,  Miss.,  January  1965. 

6.  Nutall,  C.  J.  ,  Jr.,  A  Dimensionless  Consolidation  of  WES  Data  on  the 

Performance  of  Sand  Under  Tire  Loads,  Technical  Report  No.  3-130, 
U.  S.  Army  Engineer  Waterways  Experiment  Station,  CE,  Vicksburg, 
Miss.,  December  1965. 

7.  Bekker,  M.G.  and  Janosi,  Z.  ,  Analysis  of  Towed  Pneumatic  Tires 

Moving  on  Soft  Ground,  Report  No.  RR-6,  LL-62,  U.  S.  Army  Land 
Locomotion  Laboratory,  Detroit,  Mich.,  March  I960. 

8.  Ladd,  D.  and  Ulery,  H.  ,  Jr.  ,  et  al,  Aircraft  Ground  Flotation  Investi¬ 

gation,  Parts  II  to  XIX,  AFFDL-TDR -66 -43,  Air  Force  Flight 
Dynamics  Laboratory,  Wright-Patter son  AFB,  Ohio,  August  1967. 

9.  McCrae,  J.  L.  ,  Powell,  C.  G.  ,  and  Wismer,  R.D.  ,  Performance  of 

Soils  Under  Tire  Loads,  Technical  Report  No.  3-666,  U.  S.  Army 
Engineer  Waterways  Experime nt  Station,  CE,  Vicksburg,  Miss., 
January  1965. 

10.  Smith,  M.  E.  and  Frietag,  D.  R.  ,  Deflection  of  Moving  Tires, Technical 
Report  No.  3-516,  U.  S.  Army  Engineer  Waterways  Experiment 
Station,  CE,  Vicksburg,  Miss.  ,  January  1965. 


92 


11.  Bekker,  M.  G.  ,  A  Definition  of  Soil  Trafficability,  Report  No.  41, 
Land  Locomotion  Laboratory,  Detroit,  Mich.,  Jane  1958. 


12.  Assur,  A.  ,  "Locomotion  Over  Soft  Soil  and  Snow"  Paper  No.  782F, 

Automotive  Engineering  Congress,  SAE,  Detroit,  Mich.  ,  January  1964. 

13.  Liston,  Ronald  A.  ,  et  al,  State  of  the  Art  Report  on  Off-the-Road 

Mobility  Research,  Research  Report  No.  6,  U.  S.  Army  Land  Loco¬ 
motion  Laboratory,  Detroit,  Mich.  ,  November  1966. 

14.  Frietag,  D.  R.  ,  Green,  A.  J.  ,  and  Murphy,  J.  R.  ,  Jr.  ,  Normal  Stresses 

at  the  Tire-Soil  Interface  in  Yielding  Soils,  Highway  Research  Record 
No.  74,  1964. 

15.  Wismer,  R.D.  ,  Performance  of  Soils  Under  Tire  Loads,  Technical 

Report  No.  3-666,  Report  3,  U.  S.  Army  Engineer  Waterways  Experi¬ 
ment  Station,  CE,  Vicksburg,  Miss.,  February  1966. 

16.  Aircraft  Operations  on  Unsurfaced  Soil,  Soil  Measurements  and  Analysis 

--  Project  Rough  Road  Alpha,  TR  No.  3-624,  U.  S.  Army  Engineer 
Waterways  Experiment  Station,  CE,  Vicksburg,  Miss.  ,  June  1963. 

17.  Powell,  C.  J.  and  Green,  A.  J.  ,  Performance  of  Soils  Under  Tire  Loads, 

Report  No.  3-666,  Report  2,  U.  S.  Army  Engineer  Waterways  Experi¬ 
ment  Station,  CE,  Vicksburg,  Miss.,  August  1965. 

18.  Thompson,  A.  B.  and  Smith,  M.  E.  ,  Stresses  Under  Moving  Vehicles; 

Wheeled  Vehicles,  U.  S.  Army  Engineer  Waterways  Experimentation, 
Technical  Report  No.  3-545,  i960. 

19.  Christian,  J.  T.  ,  Two-Dimensional  Analysis  of  Stress  and  Strain  in  Soils, 

Report  No.  3,  "Plane-Strain  Deformation  Analysis  of  Soil,  "  Reports 
to  DASA,  Contract  DA-22-079-eng-471,  MIT,  December  1966. 

20.  Vallabhan,  C.V.G.  and  Reese,  L.  C.  ,  Finite  Element  Method  Applied  to 

Problems  in  Stresses  and  Deformation  of  Soil,  NASA-Cr-66318,  1967. 

21.  Takagi,  S.  ,  "Plane  Plastic  Deformation  of  Soils,  "  Journal  of  Engineering 

Mechanics  Division,  ASCE,  Vol.  88,  No,  EM3,  1962. 

22.  Drucker,  D.  C.  and  Prager,  W.  ,  "S;.  .1  Mechanics  and  Plastic  Analysis 

of  Limit  Design,  "  Quarterly  Applied  Mathematics,  Vol.  10,  No.  2, 

1952.  ' .  .  . . . 


23. 


Rao,  H.A.  B.  and  Hbeg,  K.  ,  Two-Dimensional  Analysis  of  Stress  and 
Strain  in  Soils,  Report  No.  4,  "Dynamic  Response  of  Strip  Footing  on 
Elastic-Plastic  Soil,"  Reports  to  DASA,  Contract  DA-22-079-eng-471, 
MIT,  November  1967. 

24.  Jones,  R.  E.  ,  "A  Generalization  of  the  Direct  Stiffness  Method  of 

Structural  Analysis,  "  American  Institute  of  Aeronautics  and  Astro¬ 
nautics  Journal,  Vol.  2,  1964. 

25.  Ang,  A.  H.  S.  and  Rainer,  J.  H.  ,  "Model  for  Wave  Motions  in  Axi- 

Symmetric  Solids,  "  Journal  of  Engineering  Mechanics  Division,  ASCE, 
Vol.  90,  EM2,  1964. 

26.  Ang,  A.  H.  S.  ,  "Numerical  Approach  for  Wave  Motions  in  Nonlinear 

Solid  Media,  "  Matrix  Methods  in  Structural  Mechanics,  AFFDL-TR- 
66-80,  1965. 

27.  Chang,  G.  C.  and  Fair,  G.S.  ,  "Finite  Element  Approach  to  Wave 

Motions  in  Elastic-Plastic  Continua,  "  AIAA  Paper  No.  68-145,  1968. 

28.  Newmark,  N.  M.  ,  "Methods  of  Computation  for  Structural  Dynamics,  " 

Proc.  ASCE,  Vol.  85,  1959. 

29.  Sameh,  A.  H.  M.  ,  and  Ang,  A.  H.  S.  ,  Numerical  Analysis  of  Axisym- 

metric  Wave  Propagation  in  Elastic-Plastic  Layered  Media,  Civil 
Engineering  Studies,  Structural  Research  Series  No.  335,  University 
of  Illinois,  May  1968. 

30.  Kraft,  David  C.  ,  Preliminary  SingLe  Wheel  Relative  Merit  Index, 

UDRI-TR -69- 16,  University  of  Dayton  Research  Institute,  Dayton, 

Ohio,  May  1969. 


94 


appendix  I 

SINGLE  WHEEL  VERIFICATION  TESTS 
Soil  Tests  and  Preparation 
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Classification 


Both  of  the  soils  selected  for  the  test  program  have  been  used 
extensively  by  WES  in  previous  mobility  studies  and  the  classification  properties 
have  been  reported  previously.  For  comparison  purposes  the  grain  size  dis¬ 
tribution  and  limits  properties  are  given  for  the  buckshot  clay  in  Figure  54,  and 
the  grain  size  distribution  for  the  mortar  sand  is  shown  in  Figure  55. 

California  Bearing  Ratio 

The  CBR  is  a  plate  bearing  test  using  a  three-square  inch  piston  which 
is  penetrated  continuously  into  the  soil  to  a  depth  of  one-half  inch  while 
continuously  recording  the  load  resistance  with  depth.  Annular  surcharge 
weights  are  placed  around  the  piston  prior  to  its  penetration.  The  ratio  of 
the  load  at  0.  1  inch  penetration  to  that  load  supported  by  a  standard  well 
graded  crushed  gravel  multiplied  by  100  is  defined  as  the  CBR  of  the  soil. 

Cone  Penetrometer  Resistance 


The  mobility  cone  penetrometer  is  a  rod  device  having  a  thirty  degree 
cone  tip  and  has  a  cross  section  base  area  of  0.5  square  inch.  The  shaft  is 
narrowed  above  the  cone  to  minimize  the  friction  between  the  side  of  the  shaft 
and  the  hole.  The  cone  penetrometer,  which  is  pushed  into  the  soil  at  a 
standard  rate,  measures  the  resistance  to  penetration  (Cone  Index)  in  pounds 
per  square  inch.  The  Cone  Index  is  a  measure  of  soil  shear  strength  and  its 
variation  with  depth. 

Plate  Bearing  Resistance 

Plate  bearing  tests  were  conducted  on  the  soil  surface  using  three  plate 

sizes  with  diameters  of  1.4  inches,  2.8  inches,  and  4.2  inches.  The  load  in 

pounds  versus  plate  sinkage  in  inches  was  recorded  continuously  with  depth 

to  a  penetration  depth  of  approximately  ten  inches.  The  results  of  these  tests 

were  used  to  determine  the  modulus  parameters  k  ,  k  ,  and  n  in  Bekker's 

c  cp 

sinkage  equation  using  techniques  previously  detailed  (11). 

Triaxial  Tests 


Undrained  triaxial  tests  were  run  on  undisturbed  clay  specimens  taken 
from  the  test  soil  carts  for  the  soil  conditions  of  Cl  =25  and  Cl  =  50.  The  tri¬ 
axial  specimens  with  zero  confining  pressure  were  approximately  1.5  inches 
in  diameter  by  3.25  inches  high,  while  the  non-zero  confining  pressure 
specimens  were  approximately  2.8  inches  in  diameter  by  5.5  inches  high.  Each 
test  specimen  was  cut  and  trimmed  from  a  larger  test  sample.  It  was  observed 
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.  STANDARD  SIEVE  OPENING  IN  INCHES  U.S.  STANDARD  SIEVE  NUMBERS  HYDROMETER 
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that  the  undisturbed  samples  did  show  some  disturbance  as  evidenced  by 
cracks  on  the  outside  of  the  soil  sample.  The  results  of  the  triaxial  tests  on 
clay  as  shown  by  the  Mohr's  circle  relationship  are  given  in  Figure  56. 

Triaxial  tests  were  also  conducted  on  air  dried  samples  of  the  mortar 
sand.  The  specimen  size  was  2.8  inches  in  diameter  by  6.0  inches  high  and 
the  samples  were  prepared  to  a  unit  weight  of  99.  0  pcf  ±1.  0  pcf  which  was  the 
average  soil  unit  weight  in  the  soil  carts  for  the  CI=25  test  condition.  The 
results  of  the  triaxial  tests  on  the  sand  are  shown  in  Figure  57. 

Soil  Pre;  -  •  "ion_ 

The  uockshot  clay  was  processed,  placed,  and  compacted  at  predetermined 
moistures  and  densities  corresponding  to  the  desired  soil  strength.  The  soil  is 
first  passed  through  a  roller  crusher  and  then  placed  in  a  pug  mill  where  the 
soil-water  mixing  takes  place  at  a  selected  moisture  content.  The  soil  is  then 
placed  in  the  soil  carts  in  6  inch  layers  with  each  layer  compacted  by  pneumatic 
tired  rollers.  Previously  developed  empirical  relationships  between  buckshot 
clay  moisture  content  and  compactive  effort  permitted  the  soil  to  be  placed  near 
the  design  soil  conditions. 

The  mortar  sand  was  prepared  in  an  air  dried  condition.  Different  soil 
strengths  were  achieved  by  varying  the  density  of  placement  of  the  soil.  The 
sand  was  pla  ced  in  uniform  layers  which  were  screened  and  vibrated  on  the 
surface  as  the  filling  progressed.  Empirical  relationships  between  the  thickness 
of  layer,  vibratory  effect,  and  soil  density  permitted  the  sand  to  be  placed  near 
the  design  soil  conditions. 
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2  65  psi 
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NORMAL  STRESS,  psi 

Figure  56  Triaxial  Test  Results,  Buckshot  Clay 
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Triaxial  Test  Results,  \4ortar  Sand 


APPENDIX  II 

SINGLE  WHEEL  VERIFICATION  TESTS 


Soil  Uniformity  Measurements 

Mobility  cone  penetration  tests  were  conducted  at  two  meter 
intervals  along  the  soil  track  previous  to  each  test  run.  The  results 
of  these  measurements  are  given  in  Figures  58  through  70,  and  are 
identified  by  the  Test  Number  which  also  appears  in  Tables  6  and  7. 
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CONE  BIN  STATIONS,  METERS 


TEST  SECTION  CIavo=25.I 
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Figure  68  Soil  Uniformity  Measurement*,  Mortar  Sand,  Test*  10,  11,  and  12 


APPENDIX  IK 

FLOTATION  INDEX-OPERATIONS  INDEX 
COMPUTER  PROGRAM 


COMPUTER  PROGRAM  NOMENCLATURE 


A 

Tire  Contact  Area  (rigid  surface) 

AL 

Major  Axis  of  Decap 

AS 

Minor  Axis  of  Decap 

Cl 

Average  Cone  Index  of  Soil 

D 

Tire  Outside  Diameter 

DF 

Rim  Flange  Diameter 

DLTA 

Tire  Radial  Deflection 

FH 

Tire  Free  Height 

FI 

Flotation  Index 

FL 

Tire  Footprint  Length 

FW 

Tire  Footprint  Width 

LS 

Load-Strength  Ratio 

OI 

Operations  Index 

P 

Vertical  Load 

TD 

Tire  Percent  Deflection 

TP 

Tire  Ply 

TS1 

Tire  Size  Indicator 

TS2 

Tire  Size  Indicator 

TT 

Tire  Type 

W 

Tire  Width 

ZD 

Sinkage-Tire  Diameter  Ratio 

118 


FLINDX 


EFN  SOURCE  STATEMENT 


IFNI S ) 


04/23/69 


C 


9 
1C 
12 
1  A 
14 


19 

?C 

22 

24 

26 


•j  <j 

3? 

34 

35 

44 


INTEGER  TT» TP 

FLOTATION  INOEX-OPERAT I ONS  INOEX  PROGRAM-COHESIVE  SOIL 
WRITF  (  6,99) 

99  FORMAT!  1HC,14H  TT*T1 RE  TYPE/14H  TS 1-TIRE  SIZE/14H  TS2*TIRE  SIZE, 
1/13H  TP.TIRF  PLY/1 8H  O-TIRE  01 AMETE  R/15H  W«TIRE  WIDTH/ 

2? 6HT0*T !RE  DEFLECTION  PERCENT /20w  DE-FLANGE  0IAMETER/9H  P«LOAD 
3/ 1 5H  CI-CONE  INDEX/21H  0  I «OPER AT  I CNS  I NDEX/20H  FI-FLOTATION  INO 
4 =X//////////) 

J*0 
<  *2 

REAL  LS 
FORMAT  (  14) 

=ORMAT  I  I5,2F8.3»  14) 

FORMAT  (2F8.3,F6.2,F8.3) 

FORMAT  I  F 12 . 2 ) 

FORMAT  (F3.2) 

I  *0 

READ  15,9)  N 

READ  (5,10  TT,TS1,TS2,TP 
READ  (5,12)  D, W, TD,OF 
READ  (5,14)  P 
READ  (5,16)  CT 
1*1+1 

FH* ( ( D-DF ) /2 • ) 

DLTA*( (TD/100. )*FH) 

AL*(2.0)*S0RT{DLTA*(D-DLTA)  ) 

AS*(2.C)*SCRT(DLTA*(W-DLTA) ) 

F  L* ( 0 . 85 ) * ( AL ) 

IF  ( TT-3 ) 30, 30, 32 
FW=(0.93)*( AS) 

GO  TO  34 
Frf-AS 

A*(0.746)MFL)*(FW> 

L  S* ( P/ ( A*C I ) ) 

ZL*(-D.03)M0.19)*(LS) 

0  I  * ( ( ZL  )  *  5  FL ) /(P) )* 10000. 

Z0*( (01  )*(?)/( D) ) *0.0001 
F  1=  (C.C15)+n.C9)*l  ZD) 

50  FORMA- (1H0,4HTT  * ,  I  2, 10X , 5HTS1  «,F6.2, 10X,5HTS2  »,F6.2, 10X,4HTP  -, 
II  3) 

5?  FORMAT!  1H0.4H  0  *  »F6. 2*  8X» 3HP  »,F9,2, 3X.4HTD  -,F6.2, 10X,5HCI  ■  , 

lc  =>..?) 


p4  c  'V>  '-'A  -  '  1HC ,  4H  A  «,F7.2,6X,4HFL  »,F7.2> 

55  r  DR“AT ' 1  MO, 4HFI  » ,F 7.  3 , 6X , 4HUI  *,F7.3/////) 
'IRITF  (6,50  TT,TSl,TS2,TP 
(6,5?)  0,  P,  TO  ,C  I 


»•  ?  *  T  E 

•6,5'.)  a,FL 

WR'TF 

16,56)  n,o 

IF  :  •  - 

R ) 2  6 , 60 , 60 

60 

<*<  +  ? 

J  =  J  *  1 

T  ~  '  1- 

N)'.  ?,54,64 

6? 

-  r.  To 

?0 

64 

S  ro*’ 

L'T' 
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FL INDX 


EFN  SOURCE  STATEMENT 


IFNtS) 


04/23/69 


INTEGER  TT,TP 

C  FLOTATION  INDEX-OPERATIONS  INDEX  PROOR AM-COHESIONLESS  SOIL 
WRITE  (  6,99) 

99  sORMAT ( 1HO, 14H  TT-T1RE  TYPE/ 14H  TSl-TIRE  SIZE/I4H  TS2*TIRE  SIZE, 
1/13H  TP*TIRE  PLY/18H  D*TI RE  OIAM5TER/15H  W»TIRE  WIDTH/ 
226HTD*TIRE  DEFLECTION  PERCENT/20H  DF  =  FLANGE  D  I AMETER/9H  P*LOAD 
V15H  C I *CONE  INDEX/21H  OI*OPER ATIONS  INDEX/2CH  F I*FLOTAT  ION  IND 
4EX//////////) 

J*0 

<■? 

REAL  LS 

9  FORMAT  114) 

10  FORMAT  <I5,2F8.3,I4) 

12  FORMAT  (2F8.3,F6.2,F8.3) 

14  FORMAT  (F12.2) 

16  FORMAT  (F8.2) 

I«0 

19  READ  (5,9)  N 

20  READ  (5,10)  TT,TS1,TS2,TP 

22  READ  (5,12)  D, W, TO ,DF 

24  READ  (6,14)  P 

26  READ  (5,16)  CT 

I-I+l 

FH* ( ( O-OF  )/2 • ) 

DLTA-I (TO/IOO.  )*FH) 

4L«(2.0)*SORT(DLTA*(D-DLTA) ) 

AS«(2.0)*SQRT(DLTA*(W-DLTA) ) 

F  L* (0. 85)* ( AL ) 

IF  (TT-3)30, 30,32 
30  FW«(C.93)*(AS) 

GO  TO  34 
32  F w* AS 

34  A«(0.746)*(FL)*(FW) 

36  LS»(P/(A*ri)  ) 

44  7L*(-D.01)+(0,23)*(LS) 

3l«((/L)*(FL)/(P))*10000. 

ZO* ( (01  )*(P)/(0)  )*0,0001 
FI«(0.015)M3.f  9>*(Z0) 

50  FORMATl 1HC,4HTT  » , 1 2, 10X, 5HTS l  *, F6 . 2 , 10X, 5HTS2  ■ ,F6. 2 , 10X ,4HTP  », 
113) 

52  F  ORMAT ( 1HC , 4H  0  *  , F6.2, 8X, 3HP  « ,F9. 2 , 8X,4HTD  »,F6.2 , IPX ,5HCI  =  , 
1*8.2) 

54  FQt(MAT(  1H0,4H  A  *,F7.2,6X .  '*UFL  *,F7.2) 

56  FORMATl 1H0,4HF!  * ,F 7. 3 , 6X ,4H0I  *,F7.3/////) 

WRITE  (8,50)  TT,TSl,TS2,TP 
WRITE  (6,52)  0,P, TD,C I 
WRITE  (6,54)  A,FL 
WRITF  (6,56)  FI, 01 
IF  (»-'■')  26,60, 60 
60  <**+? 

J  *  J  x ' 

I  F  '  >\)*2, 64,64 

62  zn  ^  20 

64  STOP 
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APPENDIX  IV 

INCREMENTAL  STRAIN-DISPLACEMENT  RELATIONS 
FOR  STRESS  POINTS  NEAR  THE  SURFACE 


\ 


Because  there  are  no  mass  points  above  the  stress  points  that  are  on 
the  first  row  near  the  surface  of  the  soil,  the  incremental  strains  &Yrz  and 

Aez  cannot  be  calculated  by  the  usual  equations.  Approximate  expression  for 
these  incremental  strains  are  derived  below. 

Incremental  Shear  Strain 


For  the  stress  point  that  is  on  the  axis  of  symmetry,  the  equation  for 
is  the  same  as  Equation  (46),  that  is, 

Aytrz(i.j)  =  °  (A-l) 


For  the  stress  points  that  do  not  lie  on  the  axis  of  symmetry,  the 
incremental  radial  displacement,  Au,  of  a  point  on  the  soil  surface  l/2h  above 
the  stress  point  (for  example,  Point  A  in  Figure  24)  may  be  approximated  as 
the  average  of  the  incremental  radial  displacements  of  the  mass  points  located 
±h  radially  from  the  stress  point,  that  is. 


Therefore, 


(Au)A  =  -[Au(i+l,j)  +  Au(i-l,j)] 


2h 


and 


ih 


(A -2) 


(A-3 


AYrz(1,  j)  S1  C  — t^-  jj|[-Awt(U1,j)  I  (A_4) 


Incremental  Vertical  Strain 


The  vertical  stress  ofe  at  these  stress  points  are  assumed  to  be  the 
average  of  the  surrounding  stress  points;  that  is, 

a*(i,j)  =-j[atz(i+l,j-l)  +  cr‘(i+l,  j+1)  +  a‘(i-l,  j-1)  4  a‘(Ul.j+l)]  (A-5) 

and  for  the  stress  point  that  lies  on  the  axis  of  symmetry, 

a*(i,  j)  =  \  [afc(i+l,  j-1)  +  a‘(i+l,  j+1)]  (A-6) 

Z  4  Z  Z 


122 


Then  the  incremental  vertical  strain  i«  obtained  from  the  incremental 

stress-strain  relations  by  solving  for  Ae*  in  Equation  (38)  or  Equation  (43). 

z 

The  results  are  as  follows: 

If  the  stress  point  is  elastic  or  plastic  unloading, 

Ae‘(i,j)  =  o‘(i,j)-ot‘At(i,j)-Y^~[Ae‘(i,j)+Aete(i,j)l  (A-7) 

If  the  stress  point  is  plastic  loading, 

Ae'(i.j)  =|a‘ (i,j)-<r‘‘At  -  XOe*(l, j)  +  Ae'fi, j)] 

+  Qf  °»(t,j)  +  J  |  gtr(i' j)  A  S*1, jl  4  l,8(1,i)  *  Tr.(ili)AYtr,(ili>  (A‘8> 
12VF  A  2V7 

.  t  t  rt/r  Aiw)  \*" 

+  B[Acr(i,  j)  +  A€e(i,j)]j>/  X  +  2G-  qI-  +  Bj 

where  Aer,  Aeq,  Ay^.  and  ctz  are  obtained  from  Equations  (36),  (37),  (A-4), 

and  (A-5)  respectively,  and  J,  B,  Q  are  defined  by  Equations  (33),  (63),  and 

(64).  The  stresses  a  ,  aQ,  and  t  are  approximated  by  their  values  in  the 

r  o  rz 

previous  iteration. 


When  the  yield  function  f  is  significantly  above  zero  due  to  computa¬ 
tional  errors,  the  stresses  are  corrected  such  that  f  is  again  brought  near 
to  zero.  The  condition  used  for  making  the  correction  is  that  the  correction 
stress  vector  is  normal  to  the  yield  function  while  treating  each  of  the 
components  of  the  stress  tensor,  that  is,  a  >  afl»  a  »  t  »  T  ,  as  independent 
variables.  As  defined  in  Equation  (31),  r  z  rz  zr 

Yield  Function  =  f  =  al  +  /T-k 
Defining  another  yield  function  as 


F  =  J  -  (k  -  al)a 


where  F  is  a  function  of  a  ,  a  ,  t  ,  and  r  ,  and  is  greater  than  zero 

r  0  rz  zr 

when  f  is  greater  than  zero,  then  the  unit  normal  vector  to  the  function  F  is 
given  by 


(A-9) 


where  VF  is  the  gradient  of  F.  VF  has  the  following  five  components: 


— h  ftF  'if 

{7F)r  =  a^=CTr-3I  +  2a(k-aI> 
r 


*^*8  =  “e-T1  +  2a,k"ca) 

U 

-j h  1 

(VF)Z  =^"  =  az  '31  +  2a(k-aI>  (A- 10) 

z 


(VF) 


rz 


ftF  _ 

dT  rz 
rz 


(VF) 


zr 


dF 

<*T 

zr 


T 

zr 


By  defining  the  correction  vector  as 


6  a. 

i. 


j 


=  (5a  ,  6a  ,  6a  ,  6t  , 
r  h  z  rz 


6t  ) 
zr 


(A-ll) 
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I 


%  ' 


where  each  component  is  the  correction  needed  for  each  stress  component. 

In  order  to  have  the  correction  vector  6a  .  normal  to  the  yield  surface  F, 

J 

the  following  condition  is  required: 


6cr 


6a 


6 


6a 


6t 


rz 


6t 

zr 


(7F)r  (vf*)Q  (7^)z  (v£)_  (V#*). 


,  Ki' 


rz 


zr 


V*| 


From  the  definition  of  scalar  product  and  since  6a.  .  is  parallel  to  VF 


l*  J 


Therefore, 


6a..  ■  VF 

=— It — 

•J  |7F| 

6a..  6a..  •  vF 


|VF  |  |  VF  |  8 

From  Equation  (A- 10), 

|vFl3  =  [ar  -  -|l  +  2a(k-al)j3  +  [oQ  -  jl  +  2a(k-al)]8 

+  [,i-iI  +  2a(k-oa)]’+Tji  +  T»r 

After  expanding  the  squared  terms  and  simplifying, 

|  vF  | 8  =  2  J  +  12a9(k-al)a 
From  Equations  (A- 10)  and  (A- 11), 

fiCIij  ’  VF  =  ®ar  [ar  “  T1  +  "  °^]  +  6a0  [°0  "  J1  +  2a(k  -  al)  j 

+  -4I  +  2(X(k"aI)l  +  (6T-J  Trr+(fUTrr 

z  [  z  j  j  rz  rz  zr  zr 

*  6J  -  6 [(k  -  al)8]  =  6F 


(A- 12) 


(A- 13) 


(A- 14) 


(A- 15) 


(A- 16) 


(A- 17) 


6F  is  the  change  of  F  due  to  the  change  of  stresses  6q^,  and  therefore  must 
be  equal  to  the  deviation  of  F  from  the  yield  surface  (or  from  zero),  that  is 

6F  =  J  -  (k-al)8  (A- 18) 
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Substituting  Equntion.  <A-,6>.  (A-17),  (A- 18,  into  Equution  (A- 14,  giv.s 


ti=  li£l  .  %  ' ^  J-lk-oIl* 

!»*|  [tFT*~  -  IT+  Ua'tk-d)' 


From  Equations  (A- 12)  and  (A- 19), 


t  =  t  -  At  = 

ra  TZ  TZ 


K,  { 


APPENDIX  VI 


FINITE  ELEMENT  COMPUTER  PROGRAM  FOR 
ANALYTICAL  SINKAGE  ANALYSIS 


"Lumped  Parameter  Iteration  Method" 
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Some  Preliminary  Remarks  About  the  Computer  Program 

a.  The  load  curve  was  assumed  to  be  symmetrical,  that  is,  the  peak 
is  at  one-half  the  time  duration  and  the  shape  of  the  loading  portion  of  the 
curve  is  the  same  as  that  of  the  unloading  portion.  Both  portions  were 
approximated  by  20  equal  linear  segments.  The  "basic  load  curve"  in  the 
program  has  a  peak  load  of  10,  000  psf,  as  shown  in  Table  8,  and  the  curve 
is  scaled  down  proportionately  in  the  program  for  lower  peak  loads. 


Table  8 

Basic  Load  Curve 

'  i 


Dimensionless  Time 


Applied  Pressure  in  psf 


0.0 

0.0 

0.025 

-20.0 

0.050 

-110.0 

0.075 

-270.0 

0.  100 

-580.0 

0.  125 

-1050.0 

0.  150 

-1730.0 

0.  175 

-2630.0 

0.200 

-3750.0 

0.  225 

-4890. 0 

0.250 

-5890.0 

0.275 

-6770.0 

0.  300 

-7540. 0 

0.325 

-8190.0 

0.350 

-8730.0 

0.375 

-9170.0 

0.400 

-9500.0 

0.425 

-9740.0 

0.450 

-9900.0 

0.475 

-9980.0 

0.500 

-10000.0 

b.  One  of  the  input  data  items  to  the  program  is  the  time  increment,  DT. 
It  is  calculated  previous  to  use  of  the  program  using  the  stability  criterion 
discussed  in  Section  IV.  The  following  procedure  should  be  followed: 
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-  An  approximate  time  increment  ia  obtained  by  the  formula 

(At)  a  r*- — 

approx.  ZcJ 

where  h  and  Cj  are  the  grid  size  and  the  dilatational  wave  velocity,  respectively. 

-  With  the  above  approximate  time  increment  as  a  guide,  a  smaller 
time  increment,  At,  is  chosen  such  that  the  number  of  increments  in  the 

td 

entire  load  curve,  —  ,  is  a  multiple  of  40.  This  ensures  that  each  linear  seg- 
ment  of  the  load  curve  will  have  integral  number  of  time  increments. 

c.  Since  the  computer  time  required  to  run  through  a  load  curve  is  quite 
long,  the  computer  program  is  written  such  that  portions  of  the  load  curve  can 
be  run  at  separate  times  by  using  magnetic  tape  input  and  output.  This  can  be 
done  by  specifying  the  starting  load  increment  number,  LB,  and  the  ending 
load  increment  number,  LEN,  and  setting  up  the  appropriate  tapes. 

d.  The  computer  run  is  monitored  by  printing  out  the  vertical  normal 
stresses,  the  vertical  displacements,  and  the  yield  indicating  table  of  the 
region  under  the  load  in  a  format  that  has  the  alternating  pattern  similar  to 
that  of  the  finite  element  model.  The  other  stresses  and  displacements  are 
not  printed  out  because  the  volume  of  print-out  would  be  prohibitively  large. 
Also,  not  every  load  increment  is  printed  out;  the  number  of  load  increments 
that  are  skipped  is  given  by  the  Index,  ILI,  which  is  an  input  data  item. 

e.  The  numerical  value  of  each  element  of  the  yield  indicating  table 
supplies  the  following  information: 

(i)  If  -1.  0  <  YIT  <  0.  0,  the  stress  point  is  elastic. 

(ii)  If  -10.015  <  YIT  <  -10.  0,  the  yield  function  is  greater  than 

zero  but  has  not  exceeded  the  tolerance  for  yielding  (which  is  0.  015 
in  this  case),  thus  the  stress  point  is  still  considered  elastic. 

(iii)  If  0.  0  <  YIT  <  10,  000,  the  stress  point  has  yielded  and  is  loading, 
and  no  stress  correction  was  applied.  The  digits  to  the  right  of 
the  decimal  point  give  the  stress  correction  factor  (!-▼)),  which  is 
a  number  between  0.  0  and  1.  0. 

(iv)  If  30000.  0  <  YIT  <  40000.  0,  the  yield  function  has  exceeded  the 
tolerance  for  stress  correction,  and  stress  correction  has  been 
applied.  The  digits  other  than  the  ten  thousand  plate  digit  gives 
the  same  information  as  (iii). 


130 


(v)  If  20000.  0  <  YIT  <  30000. 0,  the  yield  function  is  negative  but  has 
not  gone  below  the  tolerance  (-0. 015)  for  becoming  elastic  again, 
thus  the  stress  point  is  still  considered  plastic.  The  digits  other 
than  the  ten  thousand  plate  digit  give  the  same  information  as  (ill). 

(vi)  If  40000.  0  <  YIT  <  70000. 0,  the  stress  point  is  plastic  and  un¬ 
loading.  The  digits  other  than  the  ten  thousand  plate  digit  give 
the  same  information  as  (iii). 


LIST  OF  SYMBOLS 


AL 

AP 

BE 

BE 

C 

Cl 

C2 

DFW 

DT 

DTT 

DW 

ER,  ET, 
ERZ 

FA 

FC 

FP 

G 

H 

HH 


Lame  constant  (X)  in  psf;  later  becomes  dimensionless 

(X/k) 

Soil  parameter  (a) 

A  variable  related  to  the  stress  invariants 

A  variable  in  the  plastic  relation  related  to  the  stress 
invariants  (B) 

Cohesion  (c)  in  psf 

Dilatational  wave  velocity  in  fps 

Shear  wave  velocity  in  fps 

The  value  of  the  maximum  percent  convergence  of  vertical 
displacement  between  successive  iterations  among  all  the 
mass  points 

Time  increment  (At)  in  seconds 

At 

Dimensionless  time  increment  (— — ) 

*d 

Percent  convergence  of  the  vertical  displacement  between 
successive  iterations 

EZ  Strain  increments  (Ae  ,  Ac.,  Ac  ,  and  Ay  )  at  time  t 

r  o  z  rz 

A  constant  for  controlling  whether  the  top  point  on  the  axis 
be  a  stress  point  (FA  =  -1.0)  or  a  mass  point  (FA  =  1.0) 

Yield  function 

Load  curve  adjustment  factor 

Shear  modulus  (G)  in  psf;  later  becomes  dimensionless  (G/k) 

Grid  size,  the  distance  between  a  mass  point  and  a  stress 
point  (h)  in  inches;  later  becomes  dimensionless  (h/a) 

Twice  the  grid  size  (2h) 


Index  in  the  radial  direction  (i) 


IC 

ILI 

ILP 

IT 

J 

JA 

KT 

LB 

LEN 

M 

MA 

ML 

N 

NIN 

P 

PCU 

PH 

PHI 

PIN 

PKP 

PO 


Index  for  controlling  the  particular  surrounding  stress 
point  to  be  calculated 

The  number  of  load  increments  skipped  in  the  print  out 

Index  for  controlling  which  load  increment  is  to  be 
printed  out 

Iteration  index 

Index  in  the  downward  direction  (j) 

Index  for  controlling,  during  the  final  iteration,  the  entry 
to  program  for  checking  if  the  stress  point  has  yielded 

Index  for  indicating  the  particular  linear  segment  of  the 
load  curve  that  is  being  considered 

Starting  load  increment  number  for  the  particular  computer  run 

Ending  load  increment  number  for  the  particular  computer  run 

The  number  of  grid  points  in  the  radial  direction 

Index  to  indicate  the  i  position  of  the  first  mass  point  of 
the  row 

The  i  index  for  the  last  fictitious  stress  point 

The  number  of  grid  points  in  the  downward  direction 

The  number  of  load  increments  in  a  linear  segment  of  the 
load  curve 

A  constant  in  the  plastic  relation  (Q) 

Magnitude  of  the  applied  pressure  after  each  linear  seg¬ 
ment  of  the  load  curve 

Stress  correction  factor  (t)) 

Frictional  angle  («p)  in  degrees;  later  becomes  in  radians 
Applied  pressure  increment  within  each  linear  segment 
Peak  load  of  load  pulse  in  psf 
Poisson's  ratio  (v) 
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•+ 


R 

RAT 

RH 

RHO 

RL 

SI 

SIK 

SJ 

SP 

SR,  ST.  SZ. 
SRZ 

SRS,  STS.  SZS. 
SRZS 

SRT,  STT, 

SZT,  SRZT 

SS 

TD 

TM 

TTD 

U,  UD,  UDD 

UI 

UPL 


Dimensionless  radius  (r/a) 

Stress  correction  factor  (1-tj) 

Dimensionless  radius  of  the  particular  surrounding  stress 
point  considered 

Weight  density  (p)  in  lb/ft* 

Radius  of  the  loaded  area  (a)  in  inches;  later  becomes 
in  feet 

The  current  applied  surface  pressure  that  is  prescribed 
at  the  fictitious  stress  points 

The  particular  applied  pressure  increment  being  considered 
Second  stress  invariant  of  the  stress  tensor  (J) 

Vertical  normal  stress  (o  )  x  10s;  use  for  printout  purpose 

Radial,  tangential,  vertical,  and  shear  stresses 

(o  ,  cl,  or  ,  T  )  at  time  t 
r  «  z  rz 

The  stresses  of  the  previous  iteration  at  the  surrounding 
stress  point  being  considered 

Radial,  tangential,  vertical,  and  shear  stresses 

(a  ,  0L,  O  ,  T  )  at  time  t  -  At 
'  r  fl  z  rz 

First  stress  invariant  of  the  stress  tensor  (I) 

Time  duration  of  the  load  pulse  (t^)  in  seconds 
Time  in  seconds  (t) 

Dimensionless  time  (t/t .)  at  the  end  of  each  linear 

a 

segment  of  the  load  curve 

Radial  displacement,  velocity,  and  acceleration 
(u,  u,  u  )  at  time  t 

Radial  displacement  increment  (Au)  at  time  t 

A  variable  for  indicating  unloading;  UPL  =  -1.  0  is  not 
loading  and  UPL  =  1.0  is  unloading 
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UT,  UDT, 
UDDT 

W,  WD,  WDD 

WI 

WO 

WP 

WP 

WT,  WDT, 
WDDT 

YI 

YIT 

YS 


Radial  displacement,  velocity,  and  acceleration 
(u,  u,  u  )  at  time  t-At 

Vertical  displacement,  velocity,  and  acceleration 
(w,  w,  w)  at  time  t 

Vertical  displacement  increment  (Aw)  at  time  t 
Incremental  plastic  work  done  (AW) 

Vertical  displacement  (w)  x  10*.  use  for  printout  purpose 
A  reference  plastic  work  done 

Vertical  dispiacement,  velocity,  and  acceleration 
(w,  w,  w)  at  time  t-At 

Yield  indicating  table  at  time  t 

Yield  indicating  table  at  time  t-At 

Yield  stress  in  shear  (k)  in  psf 
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M  \  _  M_  9 
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vi?  IT-:  (6,1  C  3  ) 
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WRITS  (6,100)  G ,  Pii  I 
P > -!  1  =PH  1  •  14 1  ,■  v2  7/  x  cC  « 

CC=(  ' *  —  3 1 N ( P H I  )  )*3.**0. J 
p=2.«sik;(  phi  )/cc 

1  J  1  3  V-  *  O 


yc  =  a  .  *r*rr'5  (  dm i  )  /CC 
"'SIT--  (6.1  *6)  A  p  ,  Y  3 

’••s  IT"  (6,111) 
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,\"3 !  T".  (6,107) 


137 


Cl  X3 


'•'R  i  T E  ( 6  *  i 08  )  rl  *  L‘  T  vkNi  il. 

WRITE  (6*109) 

WRITE  (6*110)  KL.TO.YS 

C  NON-O I .!c TiSl OfiAL I Z  1  <,G  ALL  P  A  iCH.-iL  T  r.  a  6  nLD  C »■» L  L UL A  [  l.  ju.  >l  i_om! I  AN  i  _  i  t  i <— 
r  '-'ILL  'JSr.'l  I N  THE  LATrlR  LOOPS 

c L=AL/YS 

r.-r- /vc 

rc,=2.*G 

H=H/RL 

MH=2.*H 

PL=RL/12. 

*'TT  =  rsT/Tr' 

C01=YS*T?>*TD*'-2.  '/  ( AH:)*FL*PL) 

003=  «P*( l.+PO)/ ( J .-t.^Po) 

DaG0/((  .  r+3.  *C03*>  P  ) 

0TG=^TT  *r.r» 

<M0  =  ~)TT<0 

r  '  >0'  ’  1  A  Y  ;-‘‘SK  LU~-'  L-jKVc.  r  Oa  Tnc  bPctir  i  ,  P.:  Lu  ••>•-*  w-^wl-I,  Lvk< 

C  I  hCR'.  "ENTS  «Kt)  ;L' ■ -0 1  I  o'^-l  i  2  i  i>  L:  I  -  i  *  •»  1 In;  a..:-i_^i\  ->r 
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.  o  t  T^  (S*1lr) 
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tr  M  J  M-M  J  W 

TTO=0. 

Fo  =  p<D/U;pr  n. 

PC=PCU( 1 ) 

DO  163  J=l*2' 

TTO=TTJ  +  O.C'2S 
JJ=J+1 

n  j  =  (  pnj  (  j  j  ) -or  j  *Pf  /fm  j  k 
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q‘J|  J)=PI/Y9 

p  r  i  1  (  j  J  \  =  A  -  • « /  v  r 

16?  VRITE  (6*116) 
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CNT°OLLIN-j  Trie.  ? f;  i  .'•! T  ..  U T  pMTTeRH* 

IF  (  r  A  •  G  T  •  •  .  ) 
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1  6  = 


REA)  (6*1^ >) 
P  t:  A  D  (  J  »  1 6  )  ) 
00  TO  1 6 S 
REA?  (r-*UO 

or  A"!  (e,i?o) 


GO  TO  1&7 
f'  T  1 


r-:Tc 

h‘T2 
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J1  =V-6 
J2  =  6-‘-’4 

I_1  =  ).?+*.’ A 
L?  =  1"7-"A 

I M  S  T  <  T  i  Wo 
•VJ1'  a;.")  '.AL-.'Ji-r-.T  - 
APPLIO  hPr.SSUK1- 

(S.1<*1) 
ILP=Li+ILI-3 
I E  =  LEN 
F  L  =  L  7  - 1 
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Ta  =  Tv/TrVn 
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<  T  =  42-.<  T 
<U  =  -1 

f I <=-P 1 N ( KT-1 ) 
r-n  TO  166 
16-  SI<  =  P1NC<T) 

1+6  '<  T ' <  T  +  < '  i 

M  =(  Tf»-Tv .)#(PC1MM!!)-PCU(<T  )  )  +PCU ( AT  )  +S I ,< 
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P 1  o  T  oLT  t  jijnL  T'.j  Zui\w/»  hiiO  Inc  Yiclj  ihuICaT Iho  .»imTi\1X  is  -»c.T  fu  ut: 
ELASTIC 
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^  y  j=i,v 
■"  =  -i 

•  "=1  .3  V'l-F  +  0.3 

.  rv>  7  T  =  1 ,  'F ,  P 

"U  T  •  J »  =n. 

VI (  I  ,  J )  =  r  • 

!  '(  I  *  J  )  =f'. 

(  I  ,  J  )  =  0  • 

!  ,  J  )  =0. 

(i  *  j  )  =c. 
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I  ,J  )  r.n, 
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nn  4  I=ML?»^1 ,7 
4  SZ(I»1)=0. 

c  STARTING  POINT  OF  THF  ITERATION  LOOP. 

JA  =  1 
I  T  =  0 

6  I  T= I T  +  l 
D  F  W  =  0  . 

F  =  F  A 

0  STARTING  POINT  FOP  THE  LOOP  I  NCR  Ex  ENT  I  MG  E«CH  RG*  GOIhj  i/U,.N««Ku 

nn  pn  j=r,nm 
F=-1.*F 

mm=?.5C01+F*o. 5 

'*F  =  M4  +  M\» 

C  STARTING  POINT  OF  THE  LOOP  FO*  CnLCUt-rt  T  I  .,o  EmCH  nG..  jTAi\TI,'<G  <>  i  t  n  IVic 
r  LIFT  pass  point 

DO  70  I  =M.M ♦  T’lF  ,  2 
TP=  1  +  1. 

IV»I-1 
R=  I  -2 
P=R#H 

r  CALCULATE,  AT  EACH  ••♦A  SO  POINT,  TriH  ALCLLERmT  I  ON-  FKU-.  The  UYgh.sIv. 
r  EQUATIONS  OF  iSuT  ION.  Trie  PART  I  CULaK  tuUnT  I  LN  T  j  Jb  c  wsPuiWj  l»iv  f  ric. 

0  LOCATION 'OF  THE  MASS  POINT. 

'JOS  =  UP(  I  ,  J  ) 

V05  =  '-'D(  I  ,J  ) 
i'S  =  U(  I  » J  ) 

PS  =  W  (  I  ,  J  ) 

I  P (  I  .GT.2)  GO  TO  x2 
UOO ( 2 , J ) =0 • 

WOO ( 2 ,J ) =C01*<  SZ <  2  » J+l >-jZ  <  2  » J-l  ) +4.*o.<Z ( 3, J )  )/rm 
GO  TO  16 

12  IF  (  J .OT . 2  •  OR •  I  .  o T  .  r;iL  )  Go  TO 

t|hh(  T  ,2  )  =0 . 

on  JO  14 

1  P  H9D(  I  ,  J  )  =C01  *  (  (  SR  (  I  P,  J  )  -  S i<  (  In,  J  )  +  SnZ  (  I  ,0+1  )-SisZ  (  I  ,  J- 1  )  )  /Hn 
1+ ( SR ( I P , J ) +SR  < I Y ♦ J ) -ST ( I P , J ) -ST ( I M . J )  ) / 2 . / R ) 

14  on!) (  1  ,J)=C01*  {  (  SZ  (  I  ,  J+l  )-SZ(  I  ,  J-l  )  +  oRA  (  IP,  J  1-jKi  (  1  .'.,wi )  )  /  tin 
1+ ( SRZ ( 1 P » J ) +SRZ ( In, J ) l/Z./R) 

CmLC’JL^T-  Tir.  0  I  SPLACi-hcriTc*  VLL>— Li  I  1  uo  Hi  I  1  l_  I  r  *\\j  1 1  I  nt 

QUADRATURE  ZGU«T I ONS . 

16  U  (  I  ,J)=*jr  (  I  » J  )  +  i9  T  l  '  o  u  T  (  I  *  J  )  +o7  T  -w^l  I*('juuT  (  i  i  »  )  /  o .  ) 

a  (  I  *J)=.vl  (  I  »  J )  +OT  F+..UT  (  I  ,  J  )+o  I  T  *(j  T  I*  (  ••u'l/i  l  1  ,^)  /  'j  .  f  • (  i  »wi)  /  c  •  ) 

U D  (  I  » J  )  =  U D T  (  I  » J )  +  l'T  T *  ( U [’■  o I  (  1  » J  )  +U.A-  (  I  ,  j  )  )  /  2  • 

vi  19  (  I  ,  J  )  =,»'  D  T  (  I  ,  J  )  +i,‘T  r*  l  Li9l>T  (  1  *  o  )  +  « •_> L  (  I  ,  -  )  )  /  2  • 

IF  (IT.'-Q.l)  GO  TO  iC 
r  AVERAGE  LIT-i  TnE  DISPLACEMENT  FkO-  THE  PRlCl.MNG  UlKaIIw, 

U(  I I  »J)+US)/2. 

V(  I  »..<)  =  (  ■'(  I  ,.l  )+'••";  >/2. 
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\'r'(  I  *  J  >  =  (  Uo(  !  »J)+L 

■,r.'(  1  ,J)  =  <•■•[;(  I  .J)-kOS)/L. 

<*  K  f.'Lr'.lL- T  “  T H1'  ■  I  •- NT  I,>1CF<  ;n-:,\iT 5 

1 0  ol  (  I  ,J)=..  (  I  .J)-cl  (  i  .J) 

V  I  (  I  » J )  -  u  (  *  » J  )  -  .  T  (  l  » J  ) 

C  D c  T L i \u  I  nc  I  OPLACt-ir  nTo  ot  Trie.  ot-Lc/i'U/  i 0  i  nc.  i-«-»  T  iw»‘<  oY  ftVLnnuiitu 
IF  (J.LT.NN)  GO  TO  19 

re  (  t  — ,  i  1  , 1  <j 

11  ‘>'1  (2»N1)=V  I(i.J)/2. 

MI ( 2  *  N1  )=0. 

r  o  T  o  n 

IF  I  (  I  - 1 » M 1  )  =  ('><!  (  I  -2  *  J  )  +’a  I  (  I  *  J )  )  /  4  • 

JI I  I  - 1  ♦  \l  )  =  (  U I  (  I  -2  » J  )  +'J  I  (  I  *  J )  )  /  4  • 

IF  (  I  <F'J.  (,-i-l  )  )  uO  TO  19 
9  ..I  (  I  + 1  >  ‘  4 1 )  =  (  :■!  I  (  l  +e  »  o  )  +  G  I  (  I  » J )  )  /  4  • 

'.i  I  (  I  + 1  *  N 1 )  =  (  U I  (  I  +2  .  J  )  +0  I  (  1  >  J )  l  /  4  . 

19  FONT  I  MU:. 

IF  IIT.LL.il  GO  Tu  i7 

t  CALF  i 1 L  ^  T  ~  T:F  Pr  PC1-  NT  COnVE  RGeNCt  OF  Tht  VORTICAL  l)  1  SPLACt.ltNT 

r  :■  c_  f  •-  £ N  T ri  i 'i  a  HQ  Trie  P K  to cii)  I  NG  1 1 h R A T 1  ON •  LuCATt  TnL  LrvRGtiT  Pc.KCtiN  I 

-*  ft  A'  o  c  A  V*-  IP  P'.ifc  LATPR  K-FF  FRCNCi . 
r'  1  ,-■.««••/>(  !  » J  )  ) 

I  p  (  O'." .  L  -  •  .'F  V  )  i_0  TO  17 

qc  v'  =  :;>v 
1L0«ILL 
IT j=IT 
I  J-  I 

J'sJ 


C  CalCULoTa  Trit  iTRi.o jio  m]  Irtc.  wJaKc-uhD I mu  oTocuo  PuINTo.  iu  CJi.TRuuu 
"  'a’HICH  STRtSS  POINT  IF  -JtiNG  COmSl Othti} 
t  I C  =  1  IS  FOR  TH-.  ONt  I  HFtD  1  aTlL  Y  olLOh 

'  IC  =  2  IS  FOR  THE  ONr  1  wVi~l>1  aTlLY  To  Tht  Li  FT 

t  10  =  3  IS  FOR  Th:;  ONI  I  .lfAEOIATF.LY  to  T  til:  RIGHT. 

1  i  t  r  -  o 


■'  =  1 


l.  =  J+l 

OH  =  P 


r  E  T  *'  R  T  I  Nu  POINT  u  F  PRuuRMr-i  i-un  CAlCOLmT  1  nG.  oTRfc.Soto  OF  cmOH  oUhROuNDING 
r  STAFFS  POINT. 


14  I C- I C+ j 

IF  (<,-.«..!)  oo  TO  64 

C  SaVE  Tut  STi<i_ewc.3  up  Trie  Pi\cCeuimu  iTLi\oTium»  ai.0  unttiv  YltLp  1  inJ  I  Cm  T  — 
C  1  riu  OaTRIX  TO  UlTlR-  Inc  «nICn  oT moo  STaaIn  RELaT i urn  To  OOC. 

SRS  =  SRU.L) 

S  T  S  =  S  T  (  <  *  L  ) 
c  7  c  =  c  7  (  K » L  ) 
c  D  Z  S  =  9  R  7.  (  K  ♦  L  ) 
t  c  ( Y I T ( Y  , L I J  2  0.0  A  ,  « 
t  STRESS  POINT  IS  FLaSTIC. 

C  CALCULATE  THE  GTRaIin  INCREmENTo  FRuri  THt  o I oPLACtotmT  I WCr\ei*ic.N  16. 

C  WHiCFi  luUaTION  To  Jot  OcPcaDo  ON  Tht  tuCArioi^  ul-  THt  oTKtoS  KOIimI. 

20  Ic  (L.GT.2)  GO  TO  24 
IF  U.GT.2)  GO  TO  el 
SZ ( 2  » 2 )  =  ( o4 I  3  » o ) +S  1  )  /  2  • 

CT  =  0. 

=?L= 0. 
to  To  76 

21  Ip  (<.GT.,,L.nn.  IC.LT.3)  GO  TO  32 
c  7  (<,?)  =  (  I  +  c,7  (  K  —  ]  ,  A)  )/?. 

r.  o  TO  ?  A 

-Jl  ..’Z  (  K.  .  2  )  =  (  G4  (  14+ 1 . 1  )  +  iZ  (  R-i  »  i  )  +62  ( IS.+  1  »  3  )  -e^>2  (  1  »  3  ) 


I 


~ - .. 


1  —  ij,z  (  <— 2  +  )/*., 

GO  TO  2  5 

3  l.  IF  ( ts  .tu.  (  iviL+  1  )  )  G  u  T  U  j  i 

2  33  SZ  ( s  *  2  )  =  ( bZ  (  K  +  l  *  1 )  +cjii  (  s.—  i  *  1 )  i-aZ  t  u+  i  »  i  )  +  jl  l  K.- 1 »  i  )  )/4» 

26  ET=  ( U I  (  K.+1  »2  )+  H  C SC. — X  *2  )  >  /  Z  .  /Krt 

E  PZ  = ( W l  (iS+l*2)— ..I(K  — 1*2)  )  /  ru  1+  ( 'J  I  (  K  »  3  )  -  U 1  ( \—  !♦*.)/<:•“  L  1  (  nt  1 » <j  )  /  c  •  )  w4 
1./3./HH 

26  rP=(UI  C  <+l  *2  )-!•!  K-l »?)  )/riri 

C7=(SZIK*2)-SZT{<.2  )-AL*(F.R  i  -£T  )  )/  (AL+  :(.  ) 

CF=4L*(F?+ET+FZI 
on  TO  30 

22  IFU.GT.2)  GO  TO  z4 
FR  =  2.*'Jl  ( J»L)/hH 
ET  =  tK 
FPZ=D. 

GO  TO  29 

24  ET=(UI  (0-l»L)+'Jl  U+±  »L  ))/£•/> <~. 

=  R  =  ( ' J I  (<+l  ,L  )-Ul  (<-l  ♦!_)  )  /nri 

I6"  (  I  T  •  GT  •  1 .OR.K.UT  ..  iL.Oft.  I <_•  LT •  3  )  Gu  Tu  ^  8 

IF  (45S(W(K+1 »L) J.GT.O.  IE-33)  GO  TO  « 

FPZ= ( U I ( <»L+1 1 -J I ( K. L-l  )  )/Hh 
GO  TO  29 

28  ERZ  =  (  W  I  ( S.+ 1  »  L  )  ~'A‘  I  (  n  —  1  * L  )  +U  I  ( .\  »  L  +  l  )  — U I  (  n  »L~  1  )  )  /  rin 

29  ZZ=  ( /-* I  (  <.«L+1  )-.vI  (n.L-1  )  l/hn 

r  CALCULATE  THE  STRESSES  FRO.'.  THh  ELASTIC  K-.um  T I  Ono. 

EE=AL*( ER+ET+FZ) 

30  SRU.L)=5RT(,<*L)+lc:+GG*c.R 
ST(<»L)=5TT  (K,L)+Lt+GG*LT 
CZ(K*L)=SZT(K»L)  +Lt +GG*tZ 
SRZ<K.L)=SRZTU»L)+G*ERZ 
GO  TO  59 

r  35  STRESS  POINT  IS  PLASTIC. 

r  CALCULATE  Thp  STRAIN  INCRF  AiENTi  KRO-i  ThL  DISFlACc  T  InCRl-v-NTu. 

<*  WHICH  EQUATION  TO  'JSC  DEPENDS  On  THU  LOCATION  Or  iHt  -■>  TF--.' iu,  ruI-,1. 

3*  UPL=— 1  a 

IF  (L.GT.c)  GC  TO  47 
T c  K.GT.2)  GO  TO  37 
SZ  (  2  »  2  )  =  (  SZ  (  3  »  3 )  +s  I  )/2. 
p  T  =  0. 
c  RZ=0 • 

GO  TU  4 i 

37  IF  (K.-.jT.OL.OR.IC.LT. 3)  Gu  TO  40 
SZ(K»2)=(SZ(<-1 »3J+SII/2. 
r,n  JO  4] 

3d  SZ  (K»2)  =  (S2(<+l»l)  +  sZ(.s—I*i  I+uZIs+ltSI+b/.  ( \  -  i  »  j ) 
i  -S7.  (  K-2  »2  )-SZ  ( K+?  » 2  ) ) / 2 • 

GO  TO  41 

40  IF  (K.eQ.  (iML+1  )  )  uu  Tu  3a 

3V  SZ  ( K  »  c  )  =  ( 'cu.  (  K+  i  »  i  )  +oZ  ( .s— a  »  i  )  +^t.  (  n+i  » j  )  l  A-  ±  •  z  )  )  /  *r  • 

41  E  T=  (  u  I  (  6-*-i  » «£ )  +U I  (K- j.  »2  )  )/*:•/  i\n 

ERZ= ( W 1  (n+1»^)-.\I(iU-i*<:  )  l/nutl  j  1  (mJ)-OI  (.v-  i  » t )  /  c.  .  — «  l  l  ,s+i » c  )  /  <.  •  )  *  h 
1./3./HH 

43  CR=('JI  (K+l  »2  )-UI  (K-1.2)  )/hn 

C  THc  NEXT  STATEMENT  CrltiCKS  The  YIeLO  1  NO  1  CM  I  No  •-ir-'  T  .x  i  A  I  ^  u  l  I  c;  a.-,  i  n  r  i  (■ 

C  THE  STRESS  POINT  IS  UNLOADING.  IT  IS  UNLOADING  IF  YIT  IS  GREaTlin  Trmh 
C  40000. 

IF  (YIT(Z.L).GT. 40000.)  GO  TO  46 

SJ  =  Sf?FT  (  (  (  SRS-5.TS )  **2+I  ST  S-SZ  ( .<  *2  )  >**c+(SZ  (  n  .  l  ) -L  *3  )  *  **+6.  *  o.-.ZS* 
1SRZS) /6. ) 

SS  =  GRS+STS+SZ(i4.2) 

HE=C03-SS/SJ/6. 
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GJ2=2.*bJ 

c7  )=CZ(<*2  )/cj2+:it 

cpt=ep+et 

F7=(SZ(.n.Z)-SzT  (6.2  )-AL*t_KT+P*bZo*  (  (  bKo*c.K+bTo*tT  +oKZi>*trtZ 
1/7.  )/SJ2+-j't*,ZRT  )  )/(GG+.AL-P*SZe*SZj) 

GO  TO  61 

46  r7=(-;Z(K,L>)_c.ZT(<,2)-AL*<PR+!;T)  )/(AL+'oCi) 

( -a  JO  A] 

47  IF(<.0T.2)  GO  TO  49 
EP  =  2.*’Jl  (  3.L  )/nH 

FT  =  ER 

C\U=0. 

GO  TO  5" 


49  E  T  =  ( U I  (,<-i  »L)+Jl  (n+i»L)  )  /  £  •  /  Km 
F°=(Ul  <<+i  ♦  L  )  -o  I  ( K.-1  *L )  )  / rtH 

CP7  =  (  WI  (6+l.L)-  '.'I  (k-1.L)+UI  U.L  +  i  >-UI  U.L-l )  ) /rth 
F Z  =  ( "( I  ( is.  .L  +  1  )-4.I  (n.L-1)  »/nri 
'1  cE=ER+LT+cZ 
*c=al*: F 

THr  MrXT  FT  ATT  MF.MT  H-iECKF  UNLOADING  aGAIN.  Ip  UNLOADING*  GO  fO 
unloading  relation. 

IF  ( Y I  T  ( s  »  L  )  •  G  T  •  40  000.  )  Go  To  53 
CaLCULaT;-.  lnt  oTKLool-  EKG.-i  Trie  HLaSTIO  NcLATIuNb. 

•:G  =  ^<b*EP+bTS*tT+o<Zb*tl<Z/2.+bz( \.L)*ci 
IF  (L.L-,.2)  GO  TO  54 
F.,  =  uR5  +  br  b  +  .E.Zb 

SJ  =  S«!RT  (  (  (  brtS-  aTv- )  **z+(  o  T o-bZb  )**<;+(  i>Zb-iKI> )  **tT&.*^Zs*atUo )  /  b.  ) 

6F=C03-SS/5J/6. 

r-J?=?.*6J 

SI  U  .L  )  =SZ  T  (  <  * L )  +GG*EZ+AE- ( bZb/S J2+oE  )  *  (  ..Q/bj*+3c.*tc  )*P 

c4  '•.r=<','0/t-  j?+ri*-*rr  )*p 

EPU.L)  =SP  T  (  <  t  L )  +(>G*GR+/\E-  ( SRb/ S J2  + at ) **lE 
FT ( \ .L )  =CT T  ( <  .L )  +CO*E  T  +  AF-  ( S Ti /  5 J2+oE  > *;vt 
F.9Z  (  s.Ll=S.<ZT(6.U+c*ER2-oRZo*..t/bJz 

GO  TO  59 

3J  'rTi\c_bb  POI N  T  Io  Uist-sMili  ns,.  UHLunuiiiu  NlcmTIwN  i  o 


3J  '~TKc.ob  POI  N T  Io  Unt-sMili ns,.  UHLunuliiu  NlcmT  iwN  i o  Joi_o. 

»  5  i  r  (l«c.ui.Z)  ju  1  c  i>o 

SZ  (  6  »L  )  =oZ  T  (  .4  »  L  )  +  Ac*’vj0i*t^Z 
5<j  SP  (  n  * L  )  =bR  T  ( .<  »L ) +Mt+s,G*  ■  K 
CT  (  K. .  L  )  =  b  T  T  (  <  *  L )  +  At.  +GG* t. T 
F?Z ( < .L ) =SRZT ( 6  *  L )  +G*EKZ 

59  CricC<  Ir  Tn-  ‘.TR-ob  POINT  hAG  YIELDED.  Tnlb  lb  00 Nr  uNLY  Fok  Trie 
c  T  VAL  T  TrP  ATI  nN.  JA  IMOFX  CONTROLb  THE  E.mTKY. 

c-p  TC  (JA.LT.2)  GO  TC  64 
C«LGUL«T  r  TriE  Y  I  U  L  .?  rONGTlU.s 
56soR  ( s  »L ) +o  T  ( >».  L )  +o<l  ( 6  »L ) 

;;j=(br<(iML)-c  T  (  ml  )  )*icii-|  jT  )  ''“trlui.listi.)  —  oi<  I )  )*** 

SP  o(=  (  okZ  (  n  *  L  )  ) 

FJ5=SSj/6.+SPQ 
SJ  =  b0RT  (  ;-Ji) 
rr=4p*cc+c J_1  . 

C:iTC.<  iF  Trie.  YiriLO  FUNCTION  16  GNEmTcK  Tpiai^  TnE  TuLcRANLt  mduVc  Zt iw» 
AN''  •  A.sc  Tiu’  ipPROPR  i  T  t  ChANGt  IN  T  i  i  Ylril.0  INDiLAT  IitC  Tm3L.Il. 

IF  (FC.GT.O.MS)  Cu  TO  63 
I c  (YFTCY.L) .LT.p.  )  GC  TO  60 
IF  (FC.GT.-n.0l6)  f;o  TO  63 

60  GO  TO  (64.bl.64) »IC 

61  YI ( K  »L) =FC 

IP  (FC.oT  »(.’•)  Y  i  ( 6  »  l.  )  =  — 1 0.  — F  L 
'■-rt  TO  64 

t'i  STrsESb  POINT  liA3  Yit-LDrO.  UnLG<  IF  IT  Io  uNLOmu I Nb. 


'VYq/  E  v  •.  .  '** 

■  .Vi*  ‘mdMMSsm 


6  3  tVO  =  Mm»L)*cR  +  oTU»L)*tT  +  SZ(i<.»L)*::Z  +  c>KZ(^.L>*U<z./<:. 

IF  (WO.GT.WR)  '.:R=WC 

IF  (WO.LT.-O.Os* ‘M.-Nu.  IC.EO.c  )  UFL=i. 
r  CALCULATE  Trip  STRESS  CORRECTION  FACTOR. 
r-SP=(  l.-AP*SS)#*? 

PH= ( S JS-SSP ) / ( 2 . *S JS+12 . *AAP*SSP > 
pf\K=PH*( (l.+6.*AAP)*SS/3.-2.*AP) 

PAT-1. -PH 

FCJ=AINT(10Qfi0.*FC) 

IF  (FC.LT.O.)  FCJ  =  4£OOCO.-FCJ 
YI U»L)=FCJ+RaT 

IF  (UPL.GT.O.)  Yi  (  n.L)  =Y1  <  ,\*  L  ) +40000.  - 
C  STRESS  CORRECTION  IS  .JlAOt  cY  Tut.  FOLcOwiNO  3  T*T  TS  oimlY  IF  Tnc 

C  YIELD  FUNCTION  IS  GREATER  THmN  A  TOL=.R*NCt. 

IF  (FC. LI .0.9ZO)  GO  TO  64 
YI  (  <  «L  )  =YI  ( <  »L  )+V 

^P<L=S°(F.L)  - 

S  T  <  I-  =  S  T  ( Y  » L ) 

3ZKL  =  SZ (K.L ) 

S  R  (  6  » L )  =SRiCL*RAT  +  PAN 
ST  (  <  * L )  =  STKL*RAT+!->AiS 
SZ  (.<»L)=3ZKL*RAT+PAiC 
SRZ(KtL)-«MT*5RZU*L) 

C  STATEMENT  64  CONTROLS  THE  o>vlTCrtlNO  To  4,\LCUL«TE  1  no  o  I  fttoou  Ai  Inc. 

C.  THREF  SURROUNDING  STRESS  POINTS. 

64  GO  TO  (65.68.70).  IC 

65  XT  =  I  - 1 
L  =  J 

RF>  =  R-H 
r,r>  TO  18 
6.0  Y  =  I  +  1 
L  =  J 

?ri  =  R+n 
GO  TO  lb 
70  CONTI  Nut 

80  CONTINUE 

Tr  ( T  T • Lc. 2 )  GO  TO  6 

'•'RITE  (6.145)  DFH*  ILP*  I  TO,  ID.  JO 

IF  (JA.EQ.Z)  GO  TO  61 

0  ChFCK  IF  THE  CONVERGENCE  IS  GCuD  LNuUGn.  RETuRN  To  CAlCuLATu  AnuTntR 
C  ITERATION  IF  NOT  ACCURATE  ENOUGH. 

IF  (OFW.GT. 0.002. AND.IT.LT. 7)  GO  TO  6 
r  IF  ACCURATE  ENOUGH.  ADJUST  JA  INDEX  AND  CALCJLmTl  FINmL  ITERATION. 
JA  =  2 
GO  TO  6 

81  CONTINUE 

C  Tut  FOLLOWING  LOOP  SAVtS  ALL  THl  0 1  SPLaClhcNTu  Fui<  Iht 

C  CALCULATION  OF  THE  NtXT  TI.4L  I NCKcnEnT  • 

F  =  -FA 

00  05  J=?.N1 
F=-1.*F 

1 . 50  0 1  +  F  *  0  •  5 

»«  —  =  ».i2— '«■* 

00  94  J=M»»,VCfp 

IF  ( I .EQ. (M+l I )  CO  TO  94 
L  =  I  + 1 

YIT ( L  t J ) =Y I ( L  » J ) 

^RT(L,J)=SR(L.J) 

STT (L.J)=ST (L.J) 

CZT(L.J)=SZ(L.J) 

SRZT(L.J)=SRZ(L,J) 
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cP( L ,J) =  S7 ( L,J)*loo. 

°4  UT  (  I  , J)  =  U(  I  ,  J) 

'•'T(  I  *  J  )  ='.*/(  I  *  J  ) 

WP(  I  ♦  J)  =V.(  I  »J)*1.CE  06 
HOT  (  I  » J ) =UD (  I  ,J) 

WOK  I  *J)=WO(  I  ,  J) 

UDDTC I ,J)=UDD( I  * J ) 

96  WODT ( I » J ) =  WL)D { I ♦ J ) 

SRT ( 1.1  >=WR 
00  180  II=MA,31,2 

180  SP( I  I *1 ) =  S2 (  I  I  ,1)*100. 

C  IF  THE  FNO  OF  A  LINEAR  SEGMENT  OF  Thi£  LuAD  CURVE  IS  REAChED »  SIlC  16 
r  CHANGED  TO  TrIE  APPLIED  PRESSuRt  INCREMENT  uF  Tht  I\caT  LINEAR  otGMENT . 
IF  ( ILL.LT.LS)  GO  TO  96 
I.  S=LS+N I  N 

i<t=kt+k" 

SI=PCU( KT ) 

5 1 K.  =  P I N  (  K.  T  ) 

96  CONTINUE 

C  PRINT  OUT  Tnb  VtRTICAL  STRcSSLG  aNU  ul  oPLAGlMlN f  o  IN  Tnu  k^oIoim  UNDlk 
C  Trie.  LO/-»L)uU  ARtA  FuR  i-ioN I  TORINO  Tut  CUMPuTiiR  RUN* 

IF  ( ILL.LT. ILP)  Gu  TO  97 

181  WRITF  (6.126)  ILL.IT.TM 
DO  183  J= 1 *NN . 2 

JJ=J+1 

WRITE  (6.FMT1)  J.(SP(I»  J ) . I = J1 . 13 .2 ) 

183  WRITE  (6.FMT2)  JJ, ( SP  (1  . JJ ) » I = J2 , 13 ,2 ) 

182  WRITE  (6,126)  ILL 
DO  184  J»2.N1,2 
JJ=J+1 

WRITE  (6»Fi>lTl)  J*(wP(l«  J  )  »  I  =  Ji  »  i.3  ,2  ) 

184  WRITE  (6.FMT2)  JJ. ( WP ( I » JJ ) » I = J2 » 13 ,2 ) 

WRITE  (6,136)  ILL, IT 

DO  186  J= 1 »NN » 2 
JJ=J+1 

WRITE  ( 6 »FMT 1 )  J,(YIT(I,  J ) , += Jl 1 1 3 ,2 ) 

186  WRITE  (6.FMT2)  JJ » ( Y I T  (  I , JJ ) , I =J2 » 13 ,2 ) 

I LP= I LP+ I L I 
WRITE  (6,130)  WR 
2C7  IF  (ILL.EO.IE)  GO  TO  201 
GO  TO  97 

2C1  WRITE  (6,125)  I LL ,  I  T  ,TM* 

DO  203  J  =  1 , N M , 2 
JJ=J+1 

WRITE  ( 6  *  FMT 1 )  J , ( SP ( I ,  J ) , I = L 1 , 25 , 2 ) 

203  WRITE  (  6  ♦  FMT  2  )  JJ ,  ( SP(  I  ,  JJ )  ,  1  =L<. ,  25 ,2  ) 

202  VQITE  (6,126)  ILL 

^0  204  J=?,N1,2 
JJ=J+1 

W9ITE  ( 6 , FMT 1 )  J,(WP(I,  J ) , I =Li , 25 ,2 ) 

204  WRITE  ( 5  » F  M  T  2 )  JJ, ( WP ( I , J J ) , I =L2 , 25 ,2 ) 

WRITE  (6,136)  ILL, IT 

DO  206  J  =  1 » N  N » 2 
JJ=J+1 

WRITE  ( 6  ,  FMT 1 )  J , <  Y I T ( I ,  J )  ,  1 =L1  ,25 ,2 ) 

206  WR  IT*  ( 6 , FMT 2 )  JJ , ( Y I T ( I , JJ )  ,  1 =L2 ,25 ,2 ) 

I E= I c+50 

97  SI=SI+SIK 

c  SAVfc  ON  TAPE  mLL  TnE  STRESSES,  i)  I  SPL''Ei*ic.NTS ,  VeLuC  I  T  I  c.o  »  ANU  aCgELcR- 
C  ATI  OHS  FOP  SUBSEQUENT  COMPUTER  RUNS. 

'••'RITE  (6,130)  WR 
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oa  «..ou=  d'  )  UT.'JLH  *l»DDT  .  vT  »i- OT«  •3['T  •i.!I  »  I 
(6*226)  <T  »L  r;  ♦  L°  *  -• !  *  c.  T  C 

99  STOP 

100  FORMAT  (  F8  •  l*F6»^*bFd.i ) 

101  FORMAT  (2F10.3»F1C.c ) 

102  FORMAT  (F8.3,F12.c<:Ib> 

103  FORMAT (lHlvlNXdS-nSUlL  PROPcAT IL3 ) 

104  FO”v  ’AT  (  22X  *  7hDl  MS  i  1  Y  *  1  /X«5r1.\nij  =*Fa0.a»1  On  c.Jo/LL--rT/c.>X».i->ii_,iiLM,A 

1  •  ■OOOLU-3 » 13X  •  3Ho  =  *  r  i.0 .  1  ♦  lCn  l^.S/Shi— FT/^i  jX»  I4mh»,  RmT  1  ^ »  x  i  a  » -n  i 

200  »,F10.2/2’X*2<?HoHFA.?  WAVS  VELOCITY  C.2  =.F10.i»7n  FT/o-C/ 

3  23X  ♦  29H0 1  L  AT  A  T I ONAL  W-VF.  VSL.  01  =*FiG.i*7n  FT/orC//> 

10b  FOR:’  -.AT  (  2  3X  *  8MCOHL  S  1  ON  .  1 5  X  *  bi  iC  =  » F 10. i » *  On  Lc- 2/c«»— r  T/t  jX  *  i.4riF  i\  I  ST  i  s 
IN  ANGLc.  »  10X  ♦  ?!  !PHI  -»F  10.1*4.-,  J_G//) 

106  F  QR;.;6  T  (  ?  3X » 2  ?HF  OR  YlcLD  v.i<  I  T-i.  1 mlP'I.-,  =  *  c  ±6  •  c /4  7X »  Siix  =*cic.o» 
1 10H  L3S/5C -CT // ) 

1 0 7  FOR’-* AT  (  IRC  » 1 9X  » 2**i 0-  >P!JTm T 1  CVi/n-  Arm  T  _R _> ) 

1C3  FOR'-’AT  ( 23X  » 1  OHoPACc  •SLbH  » 

116X»3Hri  -»F10.4*3m  I N/2  iX  » ^  ic  T  I  <-tc_  I  wC  a- f  uT  =» 

2F10. 7>4h  3c 0 / / 2 3X  »  1 1  riNUs'iwb  R  -«F  I  » 1  pX  *  j> r~. . •  •  =  »  1 4/ t_  ja  ♦  x  ih  hui'Ola  wr  _i  * 
313X»3HN  = « 1 4/ 2  .'<X  ♦  i  otiLOA!'*c.r,— FRc-  jCivjca  »  7x»4,n  ii_  =  *  I  4/ / ) 

109  FOR-’a T  ( 1UC  » 19X  *6C  *ChARACTcR  1  >T  i C  F.  A-^-Tl^  Fl.\  -.s-no  Ic.ml  I  4  I 

1\0) 

110  eoRMA  T(?**x*  1.AHLF?JC-.T  H  =  RL  =  •  r  1 0.2  ♦  3F*  I N/23X  ♦  14HT  1 1«£  =  TO  =*FiO 

1  * 2 » 4H  SEC /2 3X  « 14H‘- T  RE ES  =  <  =»F10.2*10H  Lo5/S«-FT///) 

111  FORMAT  ( IPO * 1°X  f  i  2.1LOAD  PmR«,-.c.T  Arts  ) 

113  FORMAT  (2?X.29HTK*fc  DURATION  OF  PUL  be  T  3  «*F10.^»*»n  SeC/ 

1  23X » 29HFAOI  tji  OF  Lv/*.*M.O  -'Kto  Kl  =  *F  iO.  p  ».>ri  IN/. 

2  23X » 29HPfrAx  PRESSURE  v/F  PULsc  PRP  =*ri0.i*10n  Ljo/ ou-i-  I /  /  ) 
II?  FORMAT  (1H///) 

11b  L'ORl'iAT  (6X*99h  AT  L 1  .-.cMS  1  unLeoS  c/OlL  uUi'ilhsIv/hlu 

iDI  XLNSlONti)  UIhLNMOi'.Lcs-)  l)ii'n;i\olui*L^io/  iOin  j  TIi’il. 

2  T/TD  CUKVt  Lu**0  LUr»Vc  L^n,  1  ^crtt-iicNT  l'sml. 

3  CURVE  LOAD  INCKtf’ENT// ) 

116  FORMAT  ( IX* 13*^14.9 »F 18.1 »F] 6 . 3  * f 1 6 . ^ ♦ I P 1 1° .  /»rla. 7) 

117  FORMAT  UHl*l?X*34npAPAMtTtRS  FOR  Tnli-  PaRTIOJL^R  RUN) 

UK  FORE  AT  (23X*  37FJSTART  INC  LCoO  I NCPt-;t  NT  NSolSR  L.  =  » 1 6/2 3X» :>7htNu> 
1 1  NO  LOAD  INCRE-'-NT  hUi'i-iR  LtN  -  *  I c/4 ->X»  .5  /FIsToR T  I ivO  SuRroLc.  Fkco 
2  SURE  SI  =•»  1  Pc  i  7.  7  »  i bn  (  D  l.;tUw  I  siJLc~io  ) /c^A  »  j  /nrtrtc  Ci  o  w,  rt  c  t  i\  L 

3CAu’)  INLRL.'itNT  Six  =  » lPc  1  / .  7  *  i&ri  ( Ul  i»!Li4o  I  um.1.00  )  /  *.  ja  *  1  jic  1  mi\I  1  is 
4G  T  I  M|£  »  2CX  »4HTf  i  =»Fi2.7»4h  seC//) 

11?  FORMAT  ( ?X » 23F*3TR1.SS  POINT  IS  PLASTIC//) 

120  FORMAT  (3X.iIb*5tl7.e > 

126  FORMAT  (49H1VERT ICAL  STRroS  ol  S  TRI  o-.T  lUN  (S2*iOC.)  «T  111  =»I^, 
16H»  IT  =  » 1 3 » 8H  ♦  TIMt  =.FiC.6*bn  Sr.L//) 

126  F0RMAT  I  ?f*Hl  VERTICAL  01  SPL6C.-Mr.fi  T  DI  STS!  oL’T  I  Ch  (  v:*10,.vs  lb)  AT  IL 
1L  =*!R»6H.  IT  =  » I  0  »  8H  *  T  I  »'f  s,H0.A*V»  SEC//) 

1  20  criH'fiT  (4^6) 

!30  copv&t  (?X»F14.n//) 

T’A  FORMAT  C!X*I2»2X»] 1F11.6) 

1 26  FOR.  .AT  (  40ni  Y I  i_Lw  1  hi  - 1  l—  T  i  i\o  T *■•:■  Lt  (  Y 1  T  )  n T  i  ll  —  *  i  j  »on  *  II  -  * 

1  IA//) 

:un  FORMAT  (  CF 1 1 . 2  ) 

141  FORMAT  (31  b) 

142  Fob** A  T  ( F F  •  ]  » F 1 2 •  f' ) 

14^  FORMAT  <2X.?(';17.rf*4ie  )  ) 

Fwn 

X  I  Hf.*Ap 

FX'TRv  ,UNno. 

.IJX'Z'R.  PZF  UNIT09 

JNIT^O  FI'-F  »  S(  2  )  »pF.AD Y *  INDL'T  » r  I N  * iL6  =  34  2 
FND 


*  t  p  5  J  N<  '):J 

INT^Y 

.'JNin. 

'.IV I  Tin  p  I  l 

CND 


. UNIQ. 

uni  no 

*  p  ( l )  ♦REAoY  *  I  NOU  T  *o  i  N  »  l>L.v: 


<t  ^  s  j  a 


n .  '  b  o 

.  0 

10 

(t,.-  -a.  •/, 

.'.,'  ,-.21;"  0 

3 

3  5  r. 

•: .  o 

2  O  •  ’  ; 

-llO.u 

-c  70.'. 

-  j  o  0  *  C 

-10r>"  .0 

-17  ^.C 

-t  6 2  0.0 

-p 7^0.0 

-4oyQ. 0 

-S'  6  90.0 

-0  1  7 0 . 0 

-7  phO  •  0 

-  o  7  p  u .  0 

-31 70.0 

-  -o  0  n  .  f 

-  7  7  *+  0  *  0 

-  -  3  oc  •  c 

- 9  7  0  0*0 

_  i  •>  r  'v  -  ,  r> 

x»I2.2X»AF7 

l.M 

X»12»F1  *.  •  6  • 
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